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ABSTRACT 

We have detected [Fe II] 17.94 /mi and 24.52 /im emission from a sample of 
M supergiants (// Cep, a Sco, a Ori, CE Tau, AD Per, and a Her) using the 
Texas Echelon Cross Echelle Spectrograph on NASA's Infrared Telescope Facility. 
These low opacity emission lines are resolved at R ~ 50, 000 and provide new 
diagnostics of the dynamics and thermodynamics of the stellar wind acceleration 
zone. The [Fe II] lines, from the first excited term (a 4 F), are sensitive to the 
warm plasma where energy is deposited into the extended atmosphere to form the 
chromosphere and wind outflow. These diagnostics complement previous Kuiper 
Airborne Observatory and Infrared Satellite Observatory observations which were 
sensitive to the cooler and more extended circumstellar envelopes. The turbulent 
velocities of Vt ur b — 12 — 13 km s _1 observed in the [Fe II] a 4 F forbidden lines 
are found to be a common property of our sample, and are less than that derived 
from the hotter chromospheric C II] 2325 A lines observed in a Ori, where V tur b — 
17— 19km s _1 . For the first time, we have dynamically resolved the motions of the 
dominant cool atmosph eric component d iscovered in a Ori from multi- wavelength 



radio interferometry by 



Lim et al. 



(119981 ) . Surprisingly, the emission centroids are 



quite Gaussian and at rest with respect to the M supergiants. These constraints 
combined with model calculations of the infrared emission line fluxes for a Ori 
imply that the warm material has a low outflow velocity and is located close to 
the star. We have also detected narrow [Fe i] 24.04 /mi emission that confirms 
that Fe II is the dominant ionization state in a Ori's extended atmosphere. 



Subject headings: stars: individual^ Ori, - stars: atmospheres, mass loss, winds, 
outflows - infrared: stars 
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INTRODUCTION 



M supergiants present a particular challenge in the study of mass-loss from cool 
evolved stars. For the K through mid-M spectral-types there are no working theories that 
can satisfactorily explain their observed wind properties. It has long been recognized 
that mass- loss driven by radiation pressure on dust do es not satisfy the energy-bu dget 



requirement for overcoming the gravitational potential (IHolzer fc MacGregor 



19851 ). Both 



indirect evidence from si 



( IDanchi et al 



icate dust temperatures inferred through semi-empirical modeling, 



David fe Papqularl (119901 ) . and direct evidence from infrared (IR) interferometry 



19941 ) show that the inner radius of the dominant dust features are located 



far from the stellar surface (~ 5 — 30-R*), and therefore some other mechanism is responsible 
for lifting the material out of the stellar gravitational potential. Observations reveal that 
there is insufficient hot plasma to drive thermal Parker-type winds. While mass-loss from 
some form of pulsation or convective ejection events has yet to be demonstrated, the winds 
of M supergiants often show complex structures. For example , M supergiants show multiple 



absorption in the CO 4.6 fim fundamental band (IBernat 



19811 ), and th e 12.5 urn and 20.8 



images of a Scorpii (Ml lab + B3 V) show that the dust is clumped (IMarsh et al. 



200 lh . 



/zm 



To drive the observed mass-loss rates (10~ 7 — 10~ 5 M Q yr _1 ) some process, or 
combination of processes, must substantially increase the density scale-height close to the 
star above the thermal hydrostatic value. A promising mass-loss mechanism for K and M 
stars of luminosity classes III (giants) through I (supergiants) emerg e d in the 1980's in 



the f orm of Alfven wave- driven winds (IHartmann fc MacGregor 



1980 



Hartmann fc Avrett 



19841 ) . Unlike acoustic waves and shocks which dissipate too close to the star, the long 
dissipation lengths of the non-compressive MHD waves provide a possible explanation for 
driving the observed mass-loss rates. These idealized Alfven wave-driven wind models 
(e.g., Wentzel-Kramers-Brillouin approximation) also suffer from theoretical problems that 



- 5 - 



require fine-tuni ng of the wave da mping; length to avoid terminal wind speeds in excess of 



those observed (IHolzer et al. 



19831 ). A characteristic of the 1-D Alfven wave-driven models 
was that they predicted a bloated and turbulent wind acceleration zone that was also a 
potential source of copious chromospheric emission that had been observed in many evolved 
K-M stars with the International Ultraviolet Telescope (IV E). The total Alfven energy 
fluxes and line-widths of the observed ultraviolet (UV) chromospheric emission appeared to 



be in reasonable agreement with the models if the magneti c fields were 0.1- 



Gauss), especially if area filling factors were included (e.g. 
1988k 



Hartmann et al. 



.0 m' 



1981 



(1-10 



Harper 



However, observations with spectrographs on board the Hubble Space Telescope (HST) 
revealed that this was not the case. The higher spectral resolution and higher signal-to-noise 
ratio spectra revealed that the optically thin UV emission line profiles of singly and doubly 
i onized specie s do not show the predicted trends of blue-shifted (out-flowing) centroids 



( IHarper 



200ll ). Remarkably, the low opacity line profiles of, e.g., C II] 2325 A and Si III] 



1892 A, in co ol evolved stars tend to show a small red-shift, i.e., flows down towards the 



(Carpenter et al. 


1991 


1995: 


Harper et al. 


1995) 



19951 ). For the particular case 



of the red supergiant Betelgeuse (M2 lab, a Orionis, HP 39801) multi-wavelength Very 
Large Array (VLA) radio interferometry (ILim et al.l 119981 ) revealed that the atmosphere 
is cooler and signifi cantly less ionized than the thermal structure predicted by the Alfven 



wave-drive model of 



Hartmann & Avrett 



(119841 ) . [Note that while the dominant component 



is q uite cool there i s war m /hot material embedded within it as indicated by Ha images, 



e.g. 



Hebden et a 



1987) , and HST STI S spatially resolved chromospheric spectra of C II] 



2325 A emission ( Harper fc Brownl 120061 ) . 



The wind acceleration region, in the first few radii above the photosphere, is of 
prime interest for placing empirical constraints on theories of mass-loss and is the focus 



Skinner et al. 



(|l997h . and 



Haas et al. 



of much resea r ch, e .g., I Crowley et al.l (120081) . lHarper et al.l (120051 ). iKirsch et al.l (120011 ) 



(119951 ) . This region is particularly important because 
it is where most of the energy is injected into the wind and the mechanisms responsible 
are likely to be most manifest. The energy deposited above and below the critical radius 
(~ 1.4 — 1.9i2*) controls the terminal wind speed and mass-loss rates, respectively. HST 
has revealed that the UV emission line profiles used previously are not good diagnostics of 
the wind acceleration region but are instead revealing complex chromospheric geometries 
and flows of hot plasma. This is a result of the exponential temperature sensitivity 
(oc n e exp{— he/ XkT} /^/T) of the electron collisional excitation rates for UV emission and 
the expo nential sensitivity of hydrog en ionization at chromospheric temperatures. For 



example. 



Hartmann fc Avrett 



(119841 ) show for a Ori 



~ A met + 



1 + 



6.3 x 10 3 exp{1.18 x 10 5 /T} 



n -1 



T Lya W(R/R*)VT 

where nn is the total hydrogen density (H I and H II), A met is the abundance of metal ions, 
Tiya is the H I Lya optical depth and W(R/R*) is the geometric dilution factor. These 
factors allow the total UV flux from the star to be dominated by small volumes of high 
temperature plasma. 

To study the wind acceleration in outflows we therefore seek new emission line 
diagnostics that are less sensitive to the presence of hot chromospheric material. Such 
lines naturally occur at longer wavelengths, but unfortunately the stellar photospheric 
continuum rises strongly longward of the UV and swamps potential line emission. Beyond 
the photospheric flux peak, in the mid-IR (5-25yum), the photospheric continuum has 
declined significantly and now, for M supergiants, the continuum becomes dominated 
by silicate dust emission. The mid-IR is also a good spectral region for optically thin 
emission line diagnostics. The longer wavelengths ensure much smaller Einstein decay 
rates, especially for forbidden transitions, as compared to UV and optical emission lines 
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and therefore mid-IR transitions are much less susceptible to multiple scatterings in the 
wind that would make line profile interpretation more problematic. We are interested in 
the tepid wind acceleration region so we also need to be able to distinguish its emission 
from the emission from the extended cold circumstellar envelopes (CSEs), which are known 
to emit emission lines from ground terms of atoms and singly ionized species, i.e., [O i] 
63.18 /i m and [Si III 34.81 um e mission observed with the Kuiper Airborne Observatory 



(KAO) flHaas fe Glassgold 



19931 ) a nd the [Fe III 25.99 Mm a nd [Fe II] 35.35 /im observed with 



Infrared Space Observatory (ISO) (IJusttanont et al. 



1999). Suitable candidates for wind 



acceleration diagnostics are emission lines from excited energy terms with T exc ~ 3000 K 
since the excitation energy is well in excess of the available thermal energy in the CSE 
(T gas ~ 100 — 1000 K), and are also detectable from the ground. 

In short, to study the wind acceleration in spatially unresolved spectra of M supergiants 
requires mid-IR diagnostics, a sensitive spectrograph with sufficient spectral resolution to 
resolve the line profiles, and a telescope optimized for these wavelengths at a dry site: [Fe II] 



3d a F emission, the Texas Echelon Cross Echelle Spectrograph (TEXES) (ILacy et al. 



20021 ) . and the infrared telescopes available on Mauna Kea is such a combination. 



This paper can be considered as having two main parts. The first part is centered 
around the TEXES observations of a sample of M supergiants and consists of §2 which 
describes the new [Fe II] diagnostics and their atomic data, §3 which describes the TEXES 
observations of the M supergiants, and §4 which describes the empirical properties of 
the line profiles. The second part focuses on Betelgeuse for which different independent 
observations and atmospheric models are available to help interpret the new observations. 
This part contains: §5 which discusses the details of [Fe II] 3c? 7 a 4 F line formation as well 
as that for other well studied CSE emission lines; §6 which discusses the implications of 
our findings for mass-loss mechanisms; and our conclusions which are presented in §7. Two 
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appendices are included: the first describes the procedure to flux calibrate the TEXES 
spectra, and the second describes a composite model atmosphere for Betelgeuse that is used 
to calculate mid- and far-IR line fluxes. 



2. New Infrared Diagnostics 



Figure [T] presents a partial Grotrian diagram of the two lowest energy terms of Fe II 
showing the characteristic excitation temperature defined as (Energy/fc). The [Fe II] 
25.99 fim and 35.35 /im emission lines observed with ISO are from within the ground 4s a 6 D 
term (0 — 977 cm -1 ) and probe the cool CSE. Here we use "CSE lines" to refer to emission 
from within ground energy terms, while the TEXES [Fe II] lines have a hybrid character 
being from an excited term. For some emission lines this distinction is an oversimplifcation, 
e.g., for [Fe i] 24.04 /tin where there may be a gradient in the ionization balance (§5.3.1). 
Observations of the [Fe II] 25.99 /mi and 35.35 /un kines were obtai ned with ISO-SWS at 



19991); 



spectral resolutions of R ~ 1000 for a Ori, and R ~ 250 for a Sco (lJusttanont et al. 
these resolving powers are at least a factor of 20 too low to reveal either the turbulence or 
the flow dynamics. These transitions form a ladder which ends in the ground energy leve 
(J = 9/2) and can be used to constrain the wind temperature: 35.35 fim (J^ = 5/2 — > 7/2) 
and 25.99 fim {J j{ = 7/2 -> 9/2). 



3 



An analogous ladder exists within the next term 3d 7 a A F (i.e., the first excited term: 
1872 - 3117 cm- 1 ): 24.52 /im (J^ = 5/2 -> 7/2), and 17.94 /im (J i4 = 7/2 -> 9/2). Fig. Q] 
shows these ladder s equences. These 3c? 7 a 4 F tran sitions have been observed previously by 



Kelly fe Lacy 



fll995[ ) in R ~ 10, 000 Irshell spectra (ILacv et al.lll9891 ) of the a Sco (Antares) 



1 We designate j and % as the upper and lower levels of the emission lines, respectively. 
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systeno and there is also weak coincident emission in the ISO a Ori spectrum. Since the 
emission lines present in the ISO spectra are unresolved, the emission line to continuum flux 
contrast will increase with increasing spectral resolution until the lines become resolved. 
The TEXES spectral resolution of R ~ 50, 000 provides an opportunity to detect these 
lines, resolve their line profiles at the 6 km s _1 level, and, with good signal-to-noise ratio, 
determine the emission centroid velocities to ~ 1 km s _1 . 

The 17.94 /im line lies in a spectral region with water features, both telluric and 
photospheric. A narrow telluric water line very close to the [Fe II] line can make the 
17.94 /im feature difficult to interpret, depending on the Doppler shift. In contrast, the 
24.52 /im line lies in a spectral region where the telluric attenuation varies very slowly across 
the line profile, making it more suitable for detailed emission profile analysis. In Tabled] we 
give the radiative atomic data for these diagnostics. 

There is also potential emission from between the a 4 F and a 6 D terms, namely 6.72 /im 
and 5.34 /xm which are also shown in Figure [U A characteristic of the forbidden transitions 
in the lowest terms of Fe II is that the radiative rate s within a term are stronger than the 



rates between the terms rtNussbaumer fc Storey 



19881 ) . These lines are therefore expected to 



have weaker emission than the 17.94 /xm and 24.52 /xm lines and also to sit upon a brighter, 
more complicated, stellar continuum. 



2 A subsequent discussion of these a Sco observations: "Haas, Werner, & Becklin (1996)" 
was not published (M. Haas, priv. comm.) 
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2.1. Atomic Data 



2.1.1. Radiative Data 



To utilize the high spectral resolution of the TEXES data requires accurate wavelengths, 
or wavenumbers, to establish the Doppler shifts of the line emission. We have adopted 
the most accurate laboratory wavenumbers of 407.8434 ± 0.0009 (la) cm -1 (24.52 /mi) 



and 557.5364 ± 0.0008 (la) cm 1 (17.94 /xm) which are from an ongoing project at 



Observatory to improve atomic data for forbidden iron lines (lAldenius fe Johansson 



iund 



20071 ). 



The la uncertainties correspond to 0.66 and 0.43 km s , respectively. 



Accurate Einstein decay coeffici ents (A,/ 



are a lso required if these lines are to 



provide thermodynamic constraints. iGarstand (119621 ) calculated the magnetic dipole and 
electric quadrupole transition probabilities for the 24.52 /im and 17.94 /xm lines, with the 
magnetic dipole decay probabilties completely dominating. More recent computations by 



Nussbaumer fe Storevl (119881) and t 



by 



h e IRO N Project SUPERSTRUCTURE code presented 



Quinet. Le Dourneuf fe Zeippen 



(119961 ) are both in good agreement. The latter two 



source s give AV s that are the same, which we adopt here, and these in turn are the same 



as the 



Garstangi (119621 ) values at the precision of his Table III. 



2.1.2. Collisional Data 



To establish whether the Fe II energy levels of the emitting plasma are in Local 
Thermal Equilibrium (LTE) or non-LTE , requ ires collision rates in and between the Fe II 



a 4 F and a 6 D terms. 



Pradhan fe Zhang! (119931 ) presented electron collision rate coefficients 



for forbi dden IR Fe II transi t ions t hat have an estimated uncertainty of 10-30%. Recently, 



however, iRamsbottom et al. 



)7) 



(120071 ) presented electron collision rates for temperatures that 



encompass those expected in M supergiant atmospheres that are lower by a factor of 2. 
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These uncertainties are small in comparison with estimates of hydrogen collision rates. 

Detailed collision rates for neutral hydrogen collisions have not been calcul ated, but 
estimates have been made for de-excitation rates that are of order 10~ 9 cm 3 s _1 ( lAannestad 



1973 



Bahcall fe Woli 



19681 ) . These are uncertain by an order of magnitude. If hydrogen is 
partially ionized then the total collision rates will be dominated by electron collisions, but 
in a cool photoionized stellar wind hydrogen and electron collision rates may be comparable. 
However, if the gas has a sufficiently high hydrogen density then the Fe II level populations 
have a Boltzmann (LTE) distribution and the 24.52 /zm and 17.94 /im diagnostics will then 
be insensitive to the collision rates and sensitive to the accurately known Einstein A-values. 
Large mass column densities will also tend to inhibit photon losses and drive the level 
populations towards LTE. 

In §5.2.1 we find that in the line forming region the Fe II a A F and a 6 D terms are close 
to collisional equilibrium so that current uncertainties in theoretical collisional excitation 
rates are of minor consequence to the interpretation of these mid-IR lines. 



3. TEXES OBSERVATIONS 

We have observed a sample of M supergiants, given in Table [2l with TEXES in 
high-resolution mode on the 3 m IRTF on Mauna Kea. The data described here were 
mostly obtained in 2004 October, 2005 January, and 2005 December (see Table |2] for dates). 
These observations are the longest wavelengths observed with TEXES, and were facilitated 
by a CdTe window that replaced the previous KBr window. For the long wavelength 
observations we used a 2" slit width (in the dispersion direction) to obtain the maximum 
spectral resolution and typically nodded 6" along the 17" slit to subtract the sky emission. 
The detector pixels have a linear size of 0.33", providing a fully sampled Line Spread 
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Function (LSF). The nodded star observations were interleaved with observations of a black 
thermal source and the sky. For our TEXES observations, very bright astrophysical sources 
with smooth continua suitable for flat fielding, such as asteroids, were not often available. 
We therefore used the black thermal source and sky observations to provide an approximate 
flat field. We derived a first order correction f or telluric features , and an estimate of the 



sky and telescope transmission as described in 
uncorrected for slit losses, is then given by 



Lacy et al. 



(120021 ). The source flux (F v 



F v ^ S v (obj - sky) - ^ (Tamb) (1) 
d„ (black — sky) 

where B v (T amb ) is the Planck function for the calibration source at temperature, T amb , and 
S u is the recorded signal. 

The data were collected and initially examined at the telescope using a near real-time 
reduction package written in the Interactive Data Language (IDL) which facilitated an 
efficient observing strategy. After the observing run the data were carefully optimized 
during re-calibration with the pipeline reduction software. The spatial profile on the 
detector, FWHM ~ 2.5", was used to create a template to extract the stellar spectrum. 

We observed the wavelength regions covering the [Fe II] lines discussed in §2, as well 
as some of the transitions between the a 4 F and a 6 D terms, and also the ground term 
[Fe I] 24.04 /im. Here we report on observations obtained so far; not all stars have been 
observed at all wavelengths. At long wavelengths, in high-resolution mode, only a portion 
of 4 spectral orders are recorded on the detector at one time. The recorded regions are only 
about 200 km s -1 wide and do not overlap. Because the features we observe are 20% of 
the order width, care must be taken to observe spectral features close to the center of the 
detector. An example of the spectral orders observed for the long wavelength lines in a Ori 
is shown in Figure [2j 
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3.1. Wavelength Calibration 

The wavelength scale was established by identifying telluric molecular features in 
adjacent orders. The wavelength calibration of the 24.52 /im line, however, requires special 
mention. Finding suitable telluric features near this line proved impossible, so two different 
wavelength solutions were examined. In the first case, the [Fe II] line was observed in 
2nd-order and then the filter was changed to observe a telluric feature in the 12 /im region 
in 4th-order. The grating equation was used to establish the 2nd-order wavelength solution. 
Another solution was established using telluric features about 8-orders from the emission 
feature. Both methods agree to 1km s -1 which is the level of desired accuracy for this 
science. 

3.2. Line Spread Function (LSF) 

The LSF was examined prior to mounting TEXES on the telescope using calibration 
water vapor spectra obtained from a low-pressure gas cell placed on the instrument 
entrance window. In the high resolution cross-dispersed operating mode with a 2" slit 
(in the dispersion direction) the water lines near 24/zm have emission cores that are well 
characterized by a Gaussian with R = X/ AXfwhm — 51, 700±1600. At 20/im the resolution 
is ~ 65, 000. The wings of the LSF are hard to quantify because the water lines sit upon a 
continuum. In the following analysis we adopt a Gaussian LSF with R = 52, 000 [24/A(/im)]. 

3.3. Gemini-N Observations of a Scorpii (Antares) 

During the TEXES Gemini-North engineering run in 2006 February, a spectrum was 
obtained of the [Fe II] 24.52 /im line in Antares on Feb 24. Unfortunately, no absolute 
wavelength calibration was obtained. The slit was oriented north-south, perpendicular to 
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the direction to the B 2.5 V star companion. The slit was roughly 6" long and we nodded 
the telescope east by 4" to remove sky emission and avoid potential contributions from the 
companion that lies 3" west. The spectral resolution determined by the gas cell prior to the 
observing run was R=60,000 at 18.8 /im, which suggests, via the scaling law used above, a 
spectral resolution at 24 fim of 6.4 km s _1 . 

These data were reduced in the same manner as the IRTF data. 

4. RESULTS 
4.1. Detections 

The stellar continuum near 17.94 /im is quite structured in these evolved, oxygen-rich 
M stars, and current limitations of theoretical M supergiant photospheric spectra preclude 
us from making positive detections unless the emission line stands above the adjacent 
continuum. In addition, the telluric interference near the 17.94 /im line, depending on its 
Doppler shift, can make this line difficult to analyze. We observed additional evolved K and 
M stars with widely different surface gravities and mass-loss rates to provide an empirical 
check on the structure of stellar photospheric continua. A summary of the observations and 
line detections is given in Table El At 24.5 /im the continuum is less structured but TEXES 
is not as sensitive. However, the 24.5 /im setting proved useful in establishing the presence 
of [Fe II] emission. We have dectected [Fe II] emission from all six of the M supergiants that 
we have observed. 

Figure [2] shows our first observations of a Ori at all three wavelength settings. As 
mentioned before, we record only a portion of four orders at these wavelengths. The figure 
shows how the 17.94 /mi line might be affected by a telluric water feature, while the sky 
transmission is smooth for the other two lines. The [Fe II] lines are much stronger than 
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the ground term [Fe i] line which, when the emissivities are considered, shows that iron is 
predominantly singly ionized in the extended atmosphere. 



A comparison of a Ori with CE Tau (M2 



emission from circumstellar oxygen-rich dust (ISloan fe Price 



ab) and u Cep 



(M2 



la) illustrates the effect of 



19981 ) on the line-to-continuum 



ratio. Figure [3] shows the the Fe II emission for CE Tau which is a close spectral-type proxy 
for a Ori and the line profiles are very similar. The line to continuum ratio is, however, 
substantially larger than observed in a Ori and this is, at least in part, a result of the 
much weaker (or absent) dust emission from CE Tau. /i Cep (shown in Figure HI) has much 
stronger silicate dust emission and a greatly reduced line to continuum contrast. In /i Cep 
the 17.94 jttm line is also detected but the two lines appear to have differently shaped profiles 
and the less symmetric 17.94 fim line is slightly redshifted with respect to the adopted 
V rad = +19.4 km s" 1 . This may be a result of underlying photospheric molecular features. 

We also observed AD Per (M2.5 lab), our most distant source at ~ 2kpc, and detected 
the 17.94 /im line which is shown in Figure [5j This star appears to have unusual dust 



chem i stry with carbon-ric h dust (SiC) but an oxygen-rich photosphere (ISkinner fe Whitmore 



1988 



Skinner et al. 



1990). 



The difference between a 1 Her (M5 II) and (3 Peg (M2.5 II-III) in the 17.9 /im region is 
shown in Figure [6] which reveals a narrow emission component near the stellar rest frame of 
a 1 Her. Subsequent 24.52 /im observations, not described here, confirm the [Fe II] emission. 
Both a Her's [Fe II] profiles are slightly narrower than in the more luminous counterparts. 
This figure illustrates the difficulty of identifying the emission line in this wavelength 
region, in the absence of reliable synthetic stellar spectra, when the emission is comparable 
in strength to the continuum. 

The 24.52 /im line of a Sco shown in Figure [7] is slightly wider than for the other stars 
(see Table [3]), however we find no indication of extended emission. The Antares nebulae is a 
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source of rich optical Fe II emission, e.g. 



Swings fe Preston! ( 



nearby B star companion (separation of 2.7") (IReimers et al 



9781) . which is excited by the 



2008 1. TEXES observations 



of this system may be sampling material from a slightly more extended, but still spatially 
unresolved, region than for the single stars, and this material may have different velocity 
fields. 



4.2. Properties of Line Profiles 



To characterize the observed emission, Gaussian profiles have been fit to the spectra 
and their properties are given in Table [31 a Ori was observed on several occasions and 
we have also flux calibrated these spectra as described in Appendix A. The fits to these 
individual spectra are given in Table H] which provides an indication of the reproducibility 
of the spectra and of the intrinsic variability. 

The centroid emission velocities are given with respect to the adopted stellar (center- 
of-mass) radial velocities which have typical uncertainties of at least 1 kms -1 . For M 
supergiants the center-of-mass radial velocities are not very well determined because there 



are photospheric radial velocity variations, e.g 

±3 — 5 km s _1 , that have both semi-regular (often with multiple periods o 



Jonesl (1l928l ). typically with amplitudes of 

hundreds of days 



Smith et al. 



(119891 ) . So at any given 



to several years) and short-term erratic variations, e.g. 
time the photosphere has a strong likelihood of moving with respect to the center-of-mass. 
Few M supergiants have been monitored for sufficient lengths of time to determine V ra d 
with sub-km s _1 accuracy. 

The observed emission centroid velocities are all close to the stellar center-of-mass radial 
velocities which suggests that the [Fe II] emission is not formed in the convective churning 
photospheres which undergo velocity fluctuations seen in optical absorption features. It is 
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more likely that the emitting region is larger and, or, decoupled from the variable surface 
layers. For a Ori we observe little, if any, variation with time which further supports 
this conclusion. During the 14 month period of our obs ervations a Ori's photospheric 



apparent velocity spanned a range of at least 3.7 km s 1 ( Gray 



20081 ). although we do not 



have enough epochs in common to study possible correlations. It is well known that the 



chrom osphere and wind of a Ori is decoupled from its photospheric variations ((Goldberg 



19791 ). What is apparent is that the [Fe II] emission profiles are neither blue-shifted nor 
are they flat-topped. Both of these properties exclude the possibility that the emission is 
formed with any significant outflow velocity, typically ~ 10 km s _1 for M supergiants. If the 
emission were from within such a moving flow, the absence of blue shifted profiles excludes 
the flux being formed close to the star, while the absence of a top-hat profile excludes 
the emission originating from an extended region. We conclude the emission arises from 
material with at most a small outflow velocity. Further consideration of the line formation 
in §5.1 reveals that the emission arises from close to the star. 

The observed line profiles are well resolved. Because the cores of the TEXES LSF and 
observed profiles are quite Gaussian, the intrinsic stellar most probable turbulent velocity 
{Vturb), which we assume to be isotropico , can be estimated from 



Vturb = xA&JObs) - (X km s" 1 ) 2 , (2) 



Dopp 

where X = 2.6 and 3.5 km s _1 are the corrections for the TEXES instrumental broadening 
at 17.94 and 24.52 //m, respectively. The observed [Fe II] line widths given in Tables [3] and 
H]are similar with a range of Vturb = 12 — 15 km s _1 with a Sco having the largest value. 

For a Ori, the [Fe II] line widths are similar for both the 17.94 /im and 24.52 /zm 
lines (~ 12.5 kms -1 ) and they do not change significantly between the three different 



! FWHM= 2^W2V turb = 1.6651K 



turb- 
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observing runs. The 24.04 /zm [Fe i] line is significantly narrower and indicates a different 
line forming region. Both of these forbidden line widths are significantly less than that 
derived from the chromospheric UV C II] 2325 A emission multiplet. The sky integrated 
C II] profiles observed with HST have non- Gaussia n profiles whose FWHM implies 



Vturb =19 — 21 km s 1 ( jCarpenter fe Robinson 



19971 ). Radiative transfer modeling of the 



spatially resolved HST/STIS C II] 2325 A emission reveals that these lines are slightly 
opacity broadened at the stellar limb and can be well matched with intrinsic turbulence of 
Vturh = 17 — 19 km s" 1 which changes slowly over large spatial scales: 1.5i2* < R < 3.5.R* 



jHarper fe Brown 



20061 ). This is the same spatial region over which we anticipate that the 



[Fe II] emission originates. 

The cool component of a Ori's inhomogeneous atmosphere, traced by thermal radio 
continuum observations, has now been dynamically resolved from the hot component, 
traced by UV emission lines, for the first time. The [Fe II] profiles, with their much lower 
temperature sensitivity, reflect the amplitude of the motions in the cooler plasma which 
are less than that of the hotter chromosphere. Since the cool atmospheric component 
includes the base of the wind outflow, it is these lower amplitude motions that should 
be associated with the unknown wind driving processes. For 1000-3500 K plasma these 
turbulent velocities, if interpreted as occurring on small spatial scales, imply significant 
Mach numbers. While the TEXES [Fe II] profiles are spatially unresolved (they are global 
averages), there is no evidence for outward travelling shocks moving with these velocities in 
the line forming region. 

For our TEXES [Fe II] detections the turbulent velocities are similar in all stars, 
which may not be a surprise since the s ample consists of mostly early M supergiants. 



The remarkable discovery by iLim et al. 



(119981 ) that the extended atmosphere of a Ori is 



dominated by cool, rather than hot gas as previously thought, has now been confirmed for 
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a Sco with VLA A-configuration observations made by Brown & Harper (jHarperll2009l ). 
The presence of extended cool non-chromospheric plasma with V tur b — 13 km s _1 is likely 
a common property of early M supergiants and not a rare curiosity, and deserves further 
attention. 



These [Fe II] turbulent velocities are larger than the macroturbulence requ ired to mode l 



upper 



photospheri c 12 ixm mo 



2006a) and a Ori flRvde et al 



ecular OH and H 2 absorption lines of fi Cep (IRyde et al 
2006a). T he 8 k m s _1 turbulenccl required to match a Ori's 



12 fim TEXES spectrum in IRyde et al 



model the optical: 11 kms 1 ( Gray 



12 kms 1 ( Lobel fe Dupree 



2000 



2006a) is actually small er than tha t needed to 



20081 ) and ~ 15 km s _1 



Gray 



200l| ), and near-IR: 



2000j) photospheric lines. The conclusion to be drawn from 



this is that as absorption lines are formed farther out from the star, they become less 
sensitive to the vigorous photospheric convective motions, which in turn is reflected in the 
lower macroturbulence required to match the observed line widths. At some radius where 
the extended atmosphere becomes decoupled from the photosphere, the turbulent motions 
increase once more in both the hot chromospheric and cool wind components. 



4.3. Thermal Constraints 

To place the dynamical information from the resolved line profiles in better context we 
need to establish where the emission is formed. In this subsection we will consider the most 
general formation properties and then in §5 we will consider the contribution functions of 
the TEXES and ISO CSE lines from a Ori in more detail. 

The characteristic formation temperature can be derived by assuming that the relative 
level populations of the upper (j) and lower (i) energy levels (nj,nj) can be described by a 



4 The most probable micro- and macro turbulence velocities added in quadrature. 
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Boltzmann distribution with a characteristic excitation temperature (T exc ) where 



n 



rii g% 



(3) 



The g's are the statistical weights, and Ej — Ei is the energy difference between the upper 
and lower energy levels. If the wind is isothermal and the energy levels are in thermal 
equilibrium then T exc = T gas . From the Eins tein A- values and by as suming optically 



thin emission for the ratio of the ISO fluxes, 



Justtanont et al 



(119991 ) derive excitation 



temperatures for the ground term emission of T exc ~ 1785 K and T exc ~ 1230 K for a Sco 
and a Ori, respectively. From the ratio of populations in the ground a 6 D and excited a 4 F 
terms the TEXES a Ori Fe II fluxes give T exc ~ 1520 — 1950 K. Since in a Ori the ground 
state emits at a lower characteristic temperature, this provides a lower limit for T exc for the 
a 4 F excitation region. For the ratio of fluxes within the a 4 F term we find a 3a lower limit 
of 2110 K, so the atmosphere is not isothermal. 

Figure [TO] (in Appendix B) shows the composite temperature structure for a Ori 
described in Appendix B and also shows the formation radii based on the ISO and TEXES 
temperature constraints un der these simple isothermal and optically thin assumptions. 



In the theoretical model of 



Rodgers fe Glassgoldl (119911 ) the a e D ground term emission 



originates near ~ 10-R 



lowever, we now know from the thermal radio continuum 



observations of 



Lim et al. 



(119981 ) that this temperature must occur slightly closer to the 
star, i.e., at ~ 7.4R*. The a A F emission originates interior to this at < 4.1/?* where 
the outflow velocities are expected to be small. This spatial constraint provides a partial 
explanation of why Doppler blue-shifted wind signatures are not observed. 

In summary, we find that the M supergiants share common properties in that their 
mid-IR [Fe II] line profiles appear to be quite Gaussian (rather than top-hat) and show 
no evidence of significant Doppler shifts indicative of outflow. To within the combined 
uncertainties the lines are at rest in the stellar rest frame. The line-to-continuum contrast 
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is a function of the circumstellar dust emission as expected. The characteristic excitation 
temperature places the line formation close to the star where the outflow velocities are 
expected to be low, which explains the lack of a clear wind signature. For the case of a Ori 
the line profiles do not show significant variability and the line widths are systematically 
smaller than those observed in spatially resolved UV spectra of the hotter chromosphere. 
We have dynamically resolved the turbulent motions in the dominant and pervasive cool 
atmospheric component. 



DISCUSSION: a Orionis in Context 



Where are the mid- and far-IR emission lines observed in a Ori with TEXES, 
KAO, and ISO formed? To quantify the emission contributions from different radii, a 
thermodynamic and dynamic model is required that encompasses the chromosphere, inner 
wind, and CSE. Currently no such comprehensive models exist, a Orionis provides the best 
studied example of an M supergiant and we will use the properties of this star throughout 
this section to quantify the mid- and far-IR line emission, with the reasonable assumption 
that the results will apply at some level to early M supergiants in general. 



While no comp lete atmospheric model exists, models do exist for the inner 



( IHarper et al 



2001 



HBL01 hereafter) and the outer CSE (IRodgers Glassgold 



region 



1991 



RG91 hereafter). Appendix B describes a spherical (1-D) composite dynamic (V tur b(R), 
Vwind(R)) and thermodynamic (T gas (R), p(R)) model that utilized these earlier results 
and interpolates between them. This composite model is essentially a combination of the 
spatially extended semi-emp irical model of HBL01 scaled to the recently revised stellar 



distance (IHarper et al. 



20081 ) and one of the variational thermodynamic models of RG91, 
and is referred to as the Composite Model Atmosphere. In the following calculations we use 
the cooler inner wind model which is shown as a dashed line in Fig. [10] in Appendix B. 
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The HBL01 inner region model was based upon multiwavelength spatially- resolved 
VLA data covering 0.7-6 cm combined with non-contemporaneous spatially unresolved data 
at shorter wavelengths. The HBL01 model predicts a thermal continuum flux at 100 GHz 
(0.3 cm) of 92.2 mJy which is insensitive to the wind dust emission. As part of a larger 
multiwavelength study of M supergiants and to provide a check on temporal changes in 
the extended atmosphere of a Ori, we obtained observations of a Ori, a Sco and a Her at 
100 GHz with the OVRc[] Millimeter Array and these are described next. 



5.0.1. Owens Valley Radio Observatory (OVRO) 



The OVRO observations of a Ori, a Sco, and a Her are summarized in Table El 
For a Ori, observed on 2003 November 9, four 1 GHz continuum bands were observed 
with the dual-channel analog correlator centered around 100 GHz and spanning the range 
96.5-103.5 GHz. The antennae were in the L configuration with baselines between 15-115 
m, although only 5 antennas were available during the observation. The instrumental 
gain was calibrated every 15 minutes using the quasar J0532+075. The absolute flux was 
bootstrapped from J0923+392 observations, because no planets were available during the a 
Ori transits, resulting in a 15% uncert ainty in the absolute flux scale. The calibrations were 



done with the OVRO MMA software ( 



Scoville et al 



1993) and th e images were produced 



using standard routines in Miriad (jSault. Teuben. fc Wright 



19951 ) 



For a Sco and a Her, observed in a shared track on 2004 March 30, the gains were 
calibrated using the quasars J1517-243 and J1608+104, respectively and fluxes were 
bootstrapped from these two quasars with a similar 15% uncertainty. The correlator setup 



5 The Owens Valley Radio Observatory was supported by the National Science Founda- 
tion, AST 9981546. 
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was the same as for a Ori. This track was taken in E configuration which contains several 
more extended baselines than L with baselines between 35 and 119 m and all 6 antennae 
were present throughout the track. 



by 



It is interesting to c ompare the 100 GHz fluxes with the 250 GHz fluxes measured 



Altenhoff et al. 



(119941 ) and shown in Table O At these high frequencies the earlier 
spectral- type companions of a Sco and a Her should have negligible flux contributions, e.g. 
Hjellming fc Newelll (119831 ) . The 250/100 GHz flux ratios for a Sco, a Ori, and a Her are 
4.9 ± 0.5, 4.4 ± 0.4, and 5.5 ± 0.5, respectively. When the 100 GHz fluxes are normalized 
to the product of the star's effective temperature and angular diameter squared, e.g., from 



Dyck et al. 



(119961 ). the two luminosity class lab M supergiants (a Ori and a Sco) have 
similar ratios (within 10%) while for the less luminous a Her, the ratio is about half this 
value, which may reflect it's less massive extended atmosphere. 

The HBL01 model predicted a 100 GHz flux (92 mJy) which is consis tent with the 



1982 



and references 



rather uncertain 90 GHz fluxes recorded in 1975 (INewell fc Hjellming) 
therein). Our OVRO a Ori flux was recorded about a year before the first TEXES 
observations and is slightly lower which may reflect the mean atmospheric temperature 
being slightly cooler than adopted in HBL01. 



5.1. Line Formation 

Here we examine the line formation of the forbidden excited [Fe II] and ground term 
CSE lines by computing their emission profiles and contribution functions. We assume that 
the source functions of the relevant forbidden lines are in Local Thermal Equilibrium (LTE), 
i.e., S v ~ B v (T). This can be a reasonable approximation when the particle densities in 
the line formation region are greater than the critical densities. These conditions can be 
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checked a posteriori and are discussed further in §5.2.1. 

We include the wind and turbulent velocity fields in the atomic absorption profile 
(which is assumed equal to the emission profile, i.e., complete redistribution) and compute 
the resulting spectral profiles from the formal solution of the equation of radiative transfer 
in a spherical atmosphere which runs from the upper photosphere through to the CSE. 

The adopted abundances and Einstein decay coefficients are given in Tables [1] and 
[7J For the thermal conditions in the extended envelope, stimulated emission is important 
for these mid- and far-IR transitions. The background continuous opacity is dominated 
by pure absorption and has contributions from bound-free opacity from excited levels of 
neutral species, and H and H~ free-free opacity. At these wavelengths, the bound-free 
opacity comes from hydrogenic quantum numbers n > 14, which are likely to be partially 
collisionally coupled to the continuum, so that both bound-free and H free-free opacity 
are proportional to the electron density. We performed trial n on-LTE calculation s for t he 
populations of the hydrogenic Rydberg (n, I) states, following iHummer fc Storey! (119921 ). 



and found that, at the large column densities where the continuum opacity is important for 
these mid- far-IR lines, the departure coefficients of the (n, I) levels are not significantly 
different from unity. The continuous opacity is important in the deeper layers because of 
the density sensitivity, n C ont oc n e nn, and the continuum sets the inner boundary condition 
for the line formation problem. The modelled line fluxes were obtained by measuring the 
emission above the computed local continuum and are given in Tables |5] and [HJ 
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5.2. Flux Contribution Functions 

An alternate way to estimate the line fluxes is to sum up the emission from each 
volume element in the extended atmosphere and wind, i.e., 

oo oo 

hvjiAji f n FeII rij 2 _ hi/jiAji f n FeII rij 3 



^ ^ ■-^21 / "J p^fl)^ d R = / n g " fMJ "> P esc (R)R*d\nR 

D 2 J n H n FeII D 2 J n H n FeII 

(4) 

where D is the distance to the star, hvji is the photon energy, A^ is the Einstein decay 
coefficient, nn is the total hydrogen population, n Fe ii/riH is the abundance of Fe II relative 
to hydrogen (= A Fe jj), rij/n Fe n is the ratio of the population of the upper emitting level to 
the total Fe II population, and P esc (R) is the single-flight escape probability of the photon 
emitted at radius R. This escape probability is adopted because the particle densities are 
high enough that the probability of a photon scattering and subsequently escaping is small. 
The high particle densities indicate that fine-structure level populations will be close to a 
Boltzmann distribution, so that 



n Fe ii U(T gas ) 

where U(T gas ) is the non-LTE partition function of Fe II, and Ej is the energy of the 
emitting (upper) level with respect to the ground energy level. 

The inner boundary is chosen to be deep enough that P esc (R) — > 0. The escape 
probability takes into account photons that are thermalized during line scattering and by 
continuum absorption. For low optical depths at a few stellar radii P esc is approximately 
the fraction of the sky not subtended by the star and rapidly approaches unity as the radius 
increases. Closer to the star P esc allows for photons that escape in the wings of the emission 
line when the line center optical depth is greater than unity. 

To illustrate where the emission lines originate we define a radially weighted 
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contribution function 



nH n FeII nj # 



(5) 



n H riFeii 

such that, when plotted against lni?, the area under the curve shows the relative 
contribution of different regions to the total flux. While the formal solution of the transfer 
equation provides, in principle, an exact P esc (R), here we use Eq. fl5]) to illustrate C rw . 

Convenient expressions for the escape probability have been derived for plane 
parallel geometry and cer tain spherical distributions of static scattering material 



(IKunasz fc Hummer 



1974 ). F or stellar win ds where V win d > V tur b the Sobolev escape 



probability is often employed ((Castor 



1970). Normalizing radial distances by the stellar 



radius, i.e., Z = R/R*, the continuum radial optical depth of unity occurs close to the 
stellar surface at Z cont . Here we approximate P esc (R) as the larger of the Sobolev value or 



P esc (R) = K 2 (r) E 2 ( W ) - 



1 + 



Z 



cont 



(6) 



where K 2 is the ha 



given by 



f-sky plane-parallel Doppler profile single-flight escape probability 



Hummerl (119811 ) with the mean optical depth of a static atmospherej. The E 2 



term approximates the fraction of line photons not lost to the continuum and is only 
important close to the star. The geometric term allows for stellar occultation, and the 
factor Z cont ~ 1.2 is the radius where the tangential continuum optical depth is unity. 

Figure [8] shows the normalized C rw for the TEXES and CSE lines. The narrow 
peak at Z ~ 1.5 corresponds to the high particle densities and maximum temperatures 
in the Composite Model with the sharp cut-off on the photospheric side resulting from 
continuum absorption which affects all lines in a similar fashion and eliminates any 



6 The argument of the kernel K 2 is the mean optical depth, and for a Doppler profile 
t — r oV^ where r is the static line center optical depth. 



-27- 



photospheric contributions to the emission fluxes. The decline outward of the peak is a 
com bined result of the declinin g T gas and density. Previous gene rations of theoretical, 



lartmann fe Avrett 



jl98lh and 



984), and semi-empirical models, e.g., 



Wischnewski fe Wendker 



Lobel fe Dupred (120001 ). had more extended warm chromospheres. Most of these 
can be ruled out by the observed narrowness, the absence of observed blue-shifts or wind 
broadening in the TEXES [Fe II] profiles. The small discontinuity in the figure at Z ~ 7 
reflects where the density structure of the outer wind has been merged with the inner 
density structure that is constrained by radio observations (see Appendix B). 

The TEXES and ISO [Fe II] lines clearly have different formation radii as previously 
suggested by their different characteristic excitation temperatures. The TEXES lines have 
half their emission from a region around the peak-T 5as (Z = 1.5) while the ISO lines are 
formed around Z ~ 6. The absence of wind shifted emission in the TEXES lines is a 
result of the significant contribution from the quasi-static region at chromosph eric radii 



and a bove. It is thought that this extended region, resolved with the VLA by 



Lim et al 



( 1998 ). is in the base of the wind where the velocity is small ([Harper et al.ll200ll ) and thus 



the TEXES [Fe II] lines are probing the wind, albeit at low outflow velocities. 

Fig. [8] suggests that the wind acceleration signature might be apparent in the 
ground term [Si II] 34.81 /im and [Fe II] 35.34 and 25.98 /im profiles if observed with 
sufficient spectral resolution. Note that these lines have a non-negligible flux contribution 
from within 3-R^, which, however, is often taken as the inner 



Rodgers fe Glassgold 



1991 



Haas & Glassgol 



1993 : 



the inner 


boundary condition 


Haas et al. 


1995 




Justtanont et al. 



19991 ). Clearly consideration of the material at R < 3.R* is required for detailed analysis of 
these mid- far-IR lines. 

The [O i], [C i], and [C II] lines are expected to have top-hat profiles and be centered 
close to the stellar rest frame, as the fraction of red-shifted emission occulted by the star is 
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tiny. Indeed observation of [C i] 609 /mi by iHuggins et al.l (119941 ) reveal that this line has 
similarities to mm-CO emission profiles and suggests that its very broad spatial contribution 
function also includes material traveling with the faster S2 shell (see Appendix §B2.1). 

The [Fe II] emission sits upon upper photospheric/lower chromospheric absorption and 
more reliable model fluxes require a more detailed description of the T gas structure between 
the chromospheric T sas -rise and the photosphere t han available at present. Recent VLTI 



MIDI 7.5-13.5 fiva observations (IPerrin et al 



20071 ) suggest the presence of a cool molecular 



rich region with T gas ~ 1550 K interior to 1.25-R* (corrected to the angular diameter 
adopted here). The presence of molecular material between the upper-photosphere and the 
chromospheric temperatu re peak is re miniscent of the bifurcated outer atmosphere observed 



off the solar limb in CO (lAyres 



2002) albeit on a larger fractional radial scale as befits the 



lower surface gravity of Betelgeuse. This molecular material m ay also be related to the 



Ryde et al. 



(|2002h . Future 



detection of water vapor in the outer photosphere of Arcturus by 
Atacama Large Millimeter Array (ALMA) interferometric sub-mm continuum observations 
will provide independent thermodynamic constraints down from the chromosphere towards 
the photosphere and covering this intriguing molecular region. Such ALMA data will 
complement those from the VLA that sample the chromosphere and wind. 



5.2.1. Line Source Functions 



Departures from the assumed LTE line source functions can lead to uncertainties in the 
contribution functions and line fluxes. Potential departures can be consi dered by exam ining 



the equivalent two-level atom description of the line source function Sl ( iMihalas 

J+(e' + 6) B v (T gas ) 



19781 ) 



S, 



1 + e' + rj 



(7) 
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where 

e' ee C* (1 - e~ h "' kT ) /A* (8) 

and Cji is the collisional de-excitation rate. B u (T gas ) is the Planck function at the line 
frequency, J is the mean intensity averaged over the line profile, and the other terms 9 and 
r) represent the radiative and collisional coupling bet ween the levels i,j and all other energy 



levels of the ion, i.e., many possible interactions (see 



Mihalas 



19781 for details). 



If the net rate of radiative and collisional coupling between the upper (j) [and lower 
(i)] level of the mid- and far-IR transition to all other levels, excluding i (or j), can be 
neglected, (i.e., ^ < e' and r\ <C e') then Eq. fl7|) reduces to the standard two-level atom 
description. Under these circumstances, if the downward collision rate for j — > i is higher 
than the radiative decay rate (A,j), then e' ^> 1 and the levels are in an LTE ratio with 
Sl — B{T gas ). Estimates of the critical particle densities required to establish thermal 
equilibrium between the ene rgy levels of forbidden lines are given in the compilation of 



Hollenbach fe McKed (1 198 91 ). Although the hydrogen collision rates are very uncertain, 
typically thermalization requires n# > 10 4 cm -3 , which is satisfied for radii < lOOi?*. This 
is a necessary, but not sufficient, condition for Sl — B{T gas ) because the coupling between 
other energy levels embodied by 6 and 77 can be important. 

For Fe II, collisions play a particularly important role because the first 64 fine-structure 
energy levels have the same parity and hence are coupled by collisions that compete with 
parity conserving electric quadrapole (E2) and magnetic dipole (Ml) transitions. Because 
there are so many energy levels the terms 9 and r\ are not simply evaluated, so to check the 
accuracy of the LTE source function approximation for the TEXES [Fe II] lines we have 
examined their source functions. Escape probabilities were used to approximate the net 



radiative brackets in an Fe II model with the first 769 energ y levels for the HB 



of a Ori. The atomic data are essentially those described by 



Sigut k Pradhanl (119981 ). The 



jOI model 
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ratio Sl/ B(T gas ) ~ 1 for R < 7R* with departures of at least 10% occuring in the outer line 
forming region 

Photoexcitation by chromospheric ultraviolet radiation in the allowed transition Fe II 
multiplets whose low er term is also a 4 F, e.g., Multiplet Nos. 20-31 (Aji ~ 10 5 — 10 7 s _1 ) 



rtFuhr fe Wiese 



20061 ) can lead to excitation depopulation rates in excess of the [Fe II] decay 
rates. These transitions are opaque and the depopulation rates depend on self-shielding 
which is sensitive to the wind velocity and turbulent gradients. Evaluating these rate is 
beyond the scope of the present work but we note that detailed non-LTE source functions 
are desirable for future analysis. 



5.3. Observed and Predicted Mid- and Far-IR Fluxes 

A comparison of the computed and observed fluxes for a Ori in Tables [5] & [HI reveals 
a rather unusual mismatch that is a function of formation radius. The computed TEXES 
[Fe II] fluxes are ~ 3.1 too large, the [Si II] and ISO [Fe II] emission lines are ~ 1.6 too 
large, while the lines formed at larger radii are in reasonable agreement with, or slightly 
underestimate, the more uncertain observations. 

There are different uncertainties in the calibrations and flux measurements of these 
lines (observed with TEXES, ISO, and KAO) that arise from different elements with 
their inherent uncertainties in abundance and ionization state. However, because of the 
overlapping formation radii this systematic trend is hard to explain in a simple way. These 



mid- anc 
~ 30R 



far-IR lines have large contributions inside the silicate dust shell observed at 



Danchi et al. 



19941 ) and molecular abundances and the dust /gas mass ratio are 
lower than for cooler M supergiants, suggesting that the CSE flux discrepancy is not a 
result of depletion from dust formation or molecular chemistry. The combined uncertainties 
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resulting from the observed fluxes and intrinsic variability should be < 30%, so next we 
explore other possible explanations. 



5.3.1. Ionization Balance? 
Iron 

The [Fe i] 24.04 //m emission arises from the ground term and, for a fixed ionization balance, 
it might be expected to be formed in the same region as the ISO [Fe II] 25.99 fim line. The 
observed ratio of TEXES [Fe II] to [Fe i] fluxes in a Ori sho ws that i r on is predominantly 



Rodgerd (Il990l ). This allows 



singly ionized, in agreement with theoretical calculations of 
the fluxes of the [Fe II] lines to be used as diagnostics of the amount of material in the 
extended atmosphere. The [Fe i] 24.04 /im flux can be reproduced with the Composite 
Model Atmosphere by assuming a constant Ap e i = ^O ^Apeii- However, the narrowness of 
the profile suggests that it has a stronger contribution from closer to the star where the 
turbulence is smaller and thus the ionization of iron increases with radius above the surface. 
The ionization balance in the chromosphere and inner wind region is controlled by the 
competing forces of photoionization by the strong stellar UV radiation field and radiative 
recombination. The [Fe i] 24.04 /im profile is less Gaussian and slightly asymmetric as 
compared to the TEXES [Fe II] lines and possibly has a small blue shift, in which case the 
[Fe i] may have a wind emission component. 



Castro- Carrizo et al 



(1200lh have reported a [Fe i] 24.04 /im flux from ISO grating 
spectra of 3.5 ± 0.4 x 10~ 19 W cm -2 which is significantly larger than we estimate from our 
TEXES spectra (7.7 ± xl0~ 20 W cm~ 2 ) which has a factor 40 greater spectral resolution. 
Fig. [2] reveals that there is another emission feature nearby which would be unresolved in 
ISO spectra might account, in part, for the difference in measured flux. We are unable 
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to identify this feature, but it is redshifted roughly 28 km s 1 from the peak of the 
[Fe il emiss ion 



Aoki et al 



herefore, we believe it is unlikely to be a separate component of [Fe i]. 
(119981 ) have also reported the detection of this [Fe i] line in ISO spectra in two 
carbon stars (TX PSc and WZ Cas), but it was not observed in the oxygen-rich giant 
30 Her (M6 III). These ISO observations suggest that even in this late-M giant there is 
sufficient UV flux to photoionize low ( < 13.6 eV) ionization potential metals, while in the 
carbon stars the iron is less ionized. 



from 



he predicted flux of the IS O [Fe i] 34.71/xm is consistent with the observed upper-limit 



Castro-Carrizo et al 



(j200l|). 



Other Elements 



Silicon is expected to be photoionized by the stellar UV radiation field and predominanty 
in Si II, while O I is expected to be the dominant ionization state. For carbon the 
ionization balance is more uncertain. C II dominates in the outer reaches of the CSE as the 



Galactic radiation field ionizes any remaining C I (IMamon et al. 



19881 ). Uncertainties in the 



ionization states do not appear to be the cause of the systematic discrepancies between the 
model and mid- and far-IR fluxes. 



5.3.2. Temporal Variability'; 



The fluxes given in Table [8] were observed over many years, and although there are 
hints of intrinsic variability these are at the same level as the uncertainties in the flux 
measurements. In some cases there may be off-source emission in the observing apertures 



which may mimic stellar variability (IHaas et al. 



19951 ). The TEXES observations do not 



indicate significant short time variations, and the ISO fluxes were obtained shortly after the 
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VLA observations used to construct the inner part of the atmospheric model. Therefore, 
we do not expect temporal variations sufficiently large to explain the model/observed flux 
disagreement. 

5.3.3. Temperature and Density Distribution? 

When T gas > T exc the line emissivity is rather insensitive to temperature, and for 
elements with many energy levels with the same parity as the ground state, e.g., Fe II, the 
increase in partition function further reduces the T^-sensitivity. It is only when T gas < T exc 
that the fluxes become particularly sensitive to the gas temperature. (The upper energy 
levels of the CSE lines are given in Tables [T] and [7]) . 

A combination of the assumed temperature and density distributions is the most likely 
explanation for the discrepancy between observed and model mid-IR fluxes - noting that 
the [O i] and [C II] are in reasonable agreement. The discrepancy appears to be a function 
of radius, being 3x too high for the chromosphere and wind base, 2x too high in the inner 
wind, and tending towards agreement in the outer layers. In the inner wind region the 
density structure in the Composite Model Atmosphere has been interpolated, via a simple 
wind velocity model and the equation of continuity, and is not well constrained. This could 
explain some of the flux discrepancies but not so readily the trend. 

The inner T 9a <,-structure in the Composite Model Atmosphere is derived from spatially 
resolved thermal radio emission. In the inhomogeneous atmosphere each line-of-sight 
through the stellar atmosphere intersects material of different properties: some at the 
high temperatures responsible for the UV chromospheric emission, and some much cooler 
and less ionized. The radio opacity is very sensitive to ionization (k oc n e nn, k oc n 2 e ) 
and hence has a larger contribution from the hot material than does the forbidden Fe II 
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opacity (k oc ripeii)- The over-estimation of fluxes from the excited 4 F term suggests 
that the temperature of the bulk of the plasma where the chromosphere has its largest 
filling factor is < 2500 K. Even though the radio brightness temperature inferred from the 
VLA is significantly lower than previously expected (prior to 1998) from semi-empirical 
chromospheric models, it appears that the radio brightness temperature is still greater than 
the temperature of the dominant gas component sampled by the [Fe II]. 

As one moves away from the star, the filling factor of hot chromospheric plasma 
decreases, and hence the difference between the mean temperature inferred from the VLA 
radio interferometry and the bulk gas temperature decreases. It is only by examining data 
from diagnostics with these different temperature and density dependencies that we can 
hope to unravel the complex structures in M supergiant atmospheres. We are now in an 
era where there are sufficient empirical constraints on the density, ionization, temperatures, 
and velocity fields that semi-theoretical models for the wind can be investigated. From an 
observational standpoint the largest single improvement would be to have fully resolved, 
flux calibrated, line profiles for all the CSE emission lines obtained with good pointing 
accuracy. With such profiles, both the dynamic and thermodynamic constraints of these 
important cooling channels would be realized simultaneously. 



6. Constraints on Wind Driving Mechanisms 



For these early M supergiants, radiation pressure on du st does not drive t he stellar 



1996 and the 



outflows. Most dust is located far above the stellar surface (IBester et al. 
shells are not very opaque at the wavelengths of the stellar flux peak. It has not been shown 
that radiation pressure on atoms, ions, and molecules can drive the observ ed outflows. More 



likely candidates for driving the outflo ws include some form o f pulsation (ILobel &: Dupree 



20011 ) or MHD wave propagation, e.g., 



Airapetian etal. 



(|2000h 



-35 - 



The resolved TEXES profiles provide an estimate of the energy available to drive the 
stellar wind which can be equated to that required to drive the observed mass outflow. The 
surface integrated energy flux required fro m the propagation of wave e nergy, neglecting 



wind radiative losses, can be written as (see 



F waV e (R) = AttZ 2 RI V prop CpV t 2 urb ~ M 



Holzer fc MacGregor 



19851 ) 



Z R* 2 



MV. 



2 

esc 



1 

z 



v.. 



(9) 



V esc is the surface escape speed (~ 65 km s x ) and C is a factor of or der unity tha t 



reflects line-of-sight projections and polarization of the wave motions ((Jordan 



1989 1. 



Observationally it can be inferred that the energy that drives mass-loss is mostly used to 
overcome the gravitational potential (i.e., oc V 2 SC ) with a small residual amount going into 
wind kinetic energy (i.e., oc V£,), so with Voo ~ lOkms - , the ratio (V^/Ksc) 2 is small, i.e., 
~ 0.024. 

Taking the atmospheric properties at the radius of the mid-point of the [Fe II] 
contribution functions, we have estimates for R 2 p(R), along with the measured value of 
Vturb — 12 km s _1 and M ~ 4 x 10 -6 M yr _1 . With these values the implied radial outward 
propagation velocity of the wave energy is V prop ~ 5 km s _1 at 1.5i2*. Blue-shifted emission 
from a gas outflow with this velocity is not observed in either the UV or IR emission lines. 

If the propagation speed corresponds to radial pulsations or acoustic waves then they 



are clo se to the sound speed, but they have not been observed. Note that lLobel &: Dupree 
( 20011 ) have infered non-radial pulsations with smaller amplitudes of order 1 kms~ . If 
the energy propagation is Alfvenic then the implied magnetic field fluctuations have 
SB = 0.4 Gauss. Since MHD fluctuations will damp when the amplitude approaches the 
radial field strength, B > 0.4 G. The implied plasma (3 = 8irP gas /B 2 is ~ 1, and the motions 
in the gas and magnetic field will be dynamically coupled. There are too many uncertainties 
in our current knowledge of the radial dependence of atmospheric properties of a Ori to 



be more definitive. The above arguments, namely the absence of emission indicative of 
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outward flows at 1.5i?*, suggests that either volume averaging of atmopsheric motions result 
in no outflow signature, or that the wind energy flux is carried by MHD fluctuations. The 
magnitude of the magnetic field and the order of the plasma (3 suggest that wave damping 
remains a viable mechanism to drive mass-loss in Betelgeuse. 



7. CONCLUSIONS 

We present the first resolved spectroscopy of forbidden iron emission from M 
supergiants in the 20 fim region. The TEXES spectra allow us to examine the dynamics 
and thermodynamics of the extended atm ospheres of early-type M sup ergiants. New 



accurate laboratory Ritz wavelengths from 



Aldenius &: Johansson! (120071 ). and the accurate 



and reproducible absolute wavelength scales of the TEXES spectrograph allow the [Fe II] 
17.94 /zm and 24.52 /xm emission lines to be scrutinized at the 1 kms -1 level, which is also 
the accuracy at which stellar center-of-mass radial velocities of M supergiants are known. 

Our results can be summarized as follows: 

• The [Fe II] emission lines are detected in all of our early M supergiant sample and 
the line-to-continuum flux ratios are consistent with the amount of circumstellar dust 
emission. The lines widths show little variation within our sample. 

• The a A F [Fe II] emission profiles are spectrally resolved in the TEXES spectra, and 
we have now dynamically resolved the bulk cool plasma at the base of the wind from 
the hot chromosphere. Although these lines are formed at the same radial distances 
as the hot chromosphere observed in the ultraviolet, they have smaller intrinsic line 
widths, providing clues to the atmopsheric heating and mass-loss mechanisms. 

• The emission cores of these [Fe II] lines indicate that the lines are formed close to the 
star. The absence of blue-shifted emission is in accord with low velocities expected in 
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the line forming region. 

The cool extended atmosphere has a radial velocity similar to that observed in hotter 
chromospheric UV (C II]) diagnostics at previous epochs. Neither component shows 
evidence of emission following the photospheric velocity fluctuations. 

Detailed comparison of the observed fluxes of the [Fe II] lines from a Ori and a 
composite model atmosphere are consistent with the view that Betelgeuse's extended 
atmosphere is dominated by cool gas. Early indications are that the bulk of the gas is 
even cooler than that inferred from the VLA radio interferometry, and that the filling 
factor of hot plasma declines throughout the first few stellar radii. 

We predict that spectrally resolved observations of the 25.99 /im [Fe II] line are likely 
to show a wind signature. This line is formed farther out than the a 4 F lines where the 
wind velocity is detectable with spectral resolutions of R > 50, 000, while not having 
too much contribution from the quasi-static region close to the star. This line was 
previously observed at lower spectral resolu tion with ISO, but should be observable 



with EXES (similar to TEXES) on SOFIA (IRichter et all 120061 ). 



The [Fe II] 17.94 /an and 24.52 /un line emission is co-spatial with the hot UV 
chromospheric and cool thermal radio continuum emission. The very different 
sensitivities of these diagnostics to the thermal and ionization structure are now 
beginning to constrain the filling factors of the different structural components. 

The ground term [Fe i] 24.04 /im line in a Ori is narrower than the [Fe II] which 
suggests that it is formed closer to the star where the turbulence is lower. The ratio 
of [Fe II] to [Fe i] fluxes indicates that iron is predominantly singly ionized in the 
extended atmosphere. 



In Appendix B we have constructed an extended atmosphere and wind model for 
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Betelgeuse but there are now sufficient empirical constraints to justify new theoretical 
thermodynamic and semi-empirical models that include a lower, more realistic, temperature 
boundary condition and are also constrained by the new mid- and far-IR and radio 
observations; however, this is beyond the scope of this present work. 
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Fig. 2. — TEXES observations of Betelgeuse. The stellar spectra with the flux calibration 
described in Appendix A are shown in red. An estimate of the sky transmission is shown 
in blue. The three emission lines are clearly visible and spectrally well resolved. The [Fe II] 
lines are much stronger than [Fe I] because Fe II is the dominant ionization stage in the 
extended envelope. The telluric absorption features are used to establish a wavelength scale 
accurate to < 1km s _1 . There is a telluric water line underlying the [Fe II] 17.936 /im line, 
while the shapes of the [Fe II] 24.519 and [Fe I] 24.042 /im lines are not significantly affected 
by telluric water. The underlying stellar photosphere also has shallow molecular features 
that give rise to the structured continuum. 




Fig. 3. — [Fe II] 17.94 jum and 24.52 /im emission lines in CE Tau (M2 lab). The abscissa is 
the Doppler velocity in the stellar rest frame, and the ordinate is the flux which has been 
flat-fielded, but not put onto an absolute scale. The shape of the sky transmission is shown 
beneath the spectrum. CE Tau is a spectral-type proxy for a Ori and the line centroid and 
widths for these two stars are very similar. This star has weak (or absent) dust emission. 
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Fig. 4. — [Fe II] 17.94 /im and 24.52 //m emission lines in \x Cep (M2 la). The axes are the 
same as described in Fig. [3J This supergiant has stronger silicate dust emission that a Ori 
which is a contributing factor to the low emission line to continuum ratio. 




Fig. 5.— [Fe II] 17.94 /im emission line in AD Per. The axes are the same as described in 
Fig. [31 This is the most distant M supergiant in the sample at ~ 2 kpc and has a mass-loss 
rate similar to nearby a Ori and a Sco which are ~ 200 pc. 
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Fig. 6.— Identification of [Fe II] 17.94 (Ltm emission line in a 1 Her (black) while (3 Peg (blue) 
does not show an emission feature. Subtracting the scaled (3 Peg spectrum from the a 1 Her 
that has been flat-fielded with NML Cyg reveals emission at the rest wavelength of a 1 Her. A 
Gaussian fit to this difference spectrum is shown in red and the properties are given in Table [3j 
The sky transmission for these two stars is also shown and significant additional uncertainties 
resulting from the combined telluric correction are expected. Subsequent observations of the 
24.52 fim line have confirmed this detection of [Fe II]. A color version of this figure is available 
in the electronic edition. 




Fig. 7. — [Fe II] 24.52 /im emission line in Antares. This spectrum was obtained with Gemini- 
N during an engineering run on 2006 February 24. The £1X6 S 3X6 cLS described in Fig. [31 The 
smooth profile is a Gaussian fit to the emission feature. No absolute wavelength scale was 
obtained for this spectrum. 
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Fig. 8. — Contribution functions for the emission lines observed in a Ori with TEXES, ISO 
and the KAO. The area under each curve gives the relative flux contribution. The narrow 
peak at Z ~ 1.45 is the contribution from the relatively dense inhomogeneous region that 
includes the hot chromosphere. 
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Table 1: Radiative Atomic Data for Diagnostic Infrared [Fe II] lines. 



Species 


Wavelength 3 


Wavenumber 




J low 




Jup 








w 


(cm -1 ) 


(cm -1 ) 




(cm -1 ) 




(s -1 ) b 




TEXES 


















[Fe II] 


17.9360 


557.5364 


1872.6005 


9/2 


2430.1369 


7/2 


5.84 x 10" 


-3 


[Fe II] 


24.5192 


407.8434 


2430.1369 


7/2 


2837.9803 


5/2 


3.92 x 10" 


-3 


ISO 


















[Fe II] 


25.9884 


384.7868 


0.0000 


9/2 


384.7868 


7/2 


2.13 x 10- 


-3 


[Fe II] 


35.3486 


282.8963 


384.7868 


7/2 


667.6830 


5/2 


1.57 x 10- 


-3 



"Energy levels and wavelengths (vacuum) are from I Aldenius fe Johansson! (|2007l ). 

fc Einstein A- values are from iNussbaumer fe Storevl (1988 ) and lQuinet. Le Dourneuf fe Zeippenl (1996) 



Table 2: TEXES [Fe II] and [Fe I] Observation Summary for 2004 October 5,6,11, 2005 
January 16,17, and 2005 December 9. 



Star 


Spectral- Type 


Vrad 


[Fell] 17.94 a 


[Fell] 24.52 a 


[Fel] 24.04 a 






{ 1^" yyi ci ^ 

^KHl b ) 




(//mj 




H Cep 


M2 la 


+19.4 C 


V 






a Sco 


Ml lab 


-3.5 C 




v 7 




a Ori 


M2 lab 


+20. 7 d 


v 7 


v 7 




CE Tau 


M2 lab 


+22.8 C 


v 7 


v 7 




AD Per 


M2.5 lab 


-44 c 


V 






a Her 


M5 II 


-33.1 c 


v h 




X 


/3Peg 


M2.5 II-III 


+9.1 c 


X 






Mira 


M7 III 


+63.5 C 


X 


X 




C Aur 


K4 Ib-II 


Binary 


X 






a Tau 


K5 III 


+54.3 C 


X 


X 





a Key: detection (yQ, non-detection (X), and not observed (-). 

b As judged by comparison with (3 Peg, see Fig [51 [Fe II] emission was subsequently confirmed at 24.52/im. 



tearbier-Brossat fe Fignonl |200fj). 



d Mean of 
e Mean of 



Jonesl (|1928l ) and lSanfordl (|l933h 



Barbier-Brossat fe Fignonl (|2000h and 



Mermilliod et al 



( 20081) . 
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Table 3: Line Centroid Velocities and Doppler Widths (±la a ) of the [Fe II] Emission Lines. 



Star 


Spectral- Type 


Vrad 


V b 

* cent 


V D opp (Obs) c 






(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


/i Cep 


M2 la 


+ 19.4 


1.7 + 0.2 


13.4 + 0.4 


a Sco 


Ml lab 


-3.5 


No WaveCal 


15.6 + 0.3 


a Ori 


M2 lab 


+20.7 


0.0 + 0.6 


12.5 + 0.8 


CE Tau 


M2 lab 


+22.8 


0.0 + 0.1 


12.0 + 0.2 


AD Per 


M2.5 lab 


-44 


2.0 ± 1.5 


14.4 + 1.4 


a Her 


M5 II 


-33.1 


1.7 + 0.5 


9.2±0.8 d 



a la are either the formal uncertainty of the Gaussian profile fit, or the dispersion of from multiple epoch 
measurements. 

& Centroid velocities (V cent ) are with respect to the adopted stellar center-of-mass radial velocities, V ra d- 
c Observed Doppler widths, Voopp(Obs), are defined in terms of the Full Width at Half Maximum: FWHM= 

1.665Voopp(Obs), and are uncorrected for insturmental line broadening. 
d This is a heavily blended feature, see Fig [SI and the uncertainties are dominated by systematic errors for 

this star. 
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Table 4: Properties of a Ori's TEXES [Fe II] emission lines. 



Date 


Flux a 


V C ent 


V Dopp (Obs) 


UT 


(10~ 19 WcnT 2 ) 


(km s _1 ) 


(km s _1 ) 


Fe II 24.52 fjm 








2004 Oct 05 


6.2 ±0.1 


No WaveCal 


13.2 ±0.2 


2004 Oct 06 


6.2 ±0.2 


No WaveCal 


12.4 ±0.5 


2004 Oct 11 


5.7 ±0.1 


±0.7 ±0.1 


12.8 ±0.3 


2005 Jan 16 


6.0 ±0.1 


-1.0 ±0.1 


13.3 ±0.3 


2005 Dec 09 


6.0 ±0.1 


-0.4 ±0.1 


12.9 ±0.2 


Fe II 17.94 /xm 








2004 Oct 05 


16.5 ±0.5 


0.8 ±0.2 


12.0 ±0.4 


2005 Jan 16 


16.1 ±0.1 


0.3 ±0.1 


11.7 ± 0.1 


2005 Dec 07 


16.0 ±0.2 


0.6 ±0.1 


11.4 ±0.2 


Fe I 24.04 /im 








2004 Oct 06 


0.77 ±0.03 


-2.7 ±0.2 


6.5 ±0.3 



"Flux is the emission measured above the local contiunuum and all la uncertainties are from the formal fits 
to a Gaussian profile. 
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Table 5: Properties of Betelgeuse's Infrared [Fe i] and [Fe II] lines. 



Ion 


Wavelength 


Vcent 


Vturb 


Flux 




Flux (model) 




(Vac. /im) 


(kms^ 1 ) 


(kins" 1 ) 


(W cnT 2 ) 




(W cm- 2 ) 


TEXES a 














[Fe II] 


17.9360 


+0.5 ±0.2 


12 + 0.1 


1.6 + 0.1 x 10- 


-18 


5.4 x 10~ 18 


[Fe II] 


24.5192 


+0.0 ±0.4 


13 + 0.2 


5.9 + 0.2 x 10- 


■19 


1.8 x 10~ 18 


[Fel] 


24.0423 


-2.7 + 0.2 


6.5 + 0.4 


7.7 + 0.3 x 10- 


■20 


c 7.4 x 10- 20 


ISO b 














[Fe II] 


25.9884 


unresolved 


2.8 + 0.1 x 10- 


-18 


4.4 x 10~ 18 


[Fe II] 


35.3486 


unresolved 


8.3 + 0.3 x 10- 


-19 


1.4 x 10~ 18 


[Fel] 


34.7133 


unresolved 


< 5 x 10~ 20 




c 2.2 x 10^ 20 



a TEXES fluxes are from this work. 

b ISO fluxes use the normalization described in Appendix A. 
c Assuming Ap e i = lO~ 2 Ap e . 



Table 6. OVRO 100 GHz Radio fluxes for a Sco, a Ori, and a Her. 



a Sco 



a Ori 



a Her 



Date 

Exposure Time (Hour) 

a (J2000) 

a(a) 

5 (J2000) 
a(8) 

Antenna Configuration 
Beam Size, Position Angle 
100 GHz Flux & la (mJy) 
250 GHz Flux a & la (mJy) 



2003 Nov 9 
2 

5h 55' 10.322" 

± 0.005" 
7° 24' 25.302" 
± 0.034" 
L 

8.6 x 4.3", -4.0° 15.4 x 5.0", -11.0° 4.9 x 4.1", +75.8° 
70 ± 1.8 80.4 ±3.7 19 ± 0.49 

345 ± 34 351 ± 25 104 ± 10 



2004 Mar 30 
2 

16h 29' 24.492" 
± 0.002" 

-26° 25' 54.676" 
± 0.005" 
E 



2004 Mar 30 

5 

17h 14' 38.862" 
± 0.002" 

14° 23' 25.611" 
± 0.002" 
E 



a 250 GHz fluxes are from 



Altenhoff et al 



fll994h . 
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Table 7: Atomic Data for CSE emission lines, and adopted abundances for a Ori. 
Transition Wavelength A,j a Abundance E up Source for Abundance 



(fim) (s^ 1 ) (Rel. to H) (cm- 1 ) 

[Fe I] 24.04 2.51d-03 3.0d-05 415.933 Carr et al. (2000 1 

[Fe II] 25.99 2.13d-03 3.0d-05 384.790 Carr et al. (2000 1 

[Fe II] 35.35 1.57d-03 3.0d-05 667.683 Carr et al. (2000 1 

[0 1] 63.18 8.91d-05 6.3d-04 158.265 Lambert et al. (1984 1 

[O I] 145.5 1.75d-05 6.3d-04 226.977 Lambert et al. (1984) 

[Sill] 34.81 2.13d-04 3.8d-05 287.24 Rodgers fe Glassgold (1991) 

[C I] 609.7 7.88d-08 2.5d-04 16.40 Lambert et al. (1984) 

[C II] 157.7 2.30d-06 2.5d-04 63.42 Lambert et al. (1984) 



"Einstein A-values are from NIST^ except for Fe I ([Brown &: Evenson 



19951) . 



tRalchenko, Yu., Kramida, A.E., Reader, J. and NIST ASD Team (2008). NIST Atomic Spectra Database 



(version 3.1.4), [Online]. Available: http://physics.nist.gov/asd3 [2008, April 16]. National Institute of 
Standards and Technology, Gaithersburg, MD. 



-63- 



Table 8: Other observed forbidden line fTuxes a for a Ori and computed flxues from the Com- 
posite Model Atmosphere described in Appendix B. 



Transition 


Flux (observed) 


Flux (model) 


Reference for Observed Fluxes 




(W cm" 2 ) 


(W cm" 2 ) 




[0 I] 63.18 /mi 


2.4 ±0.2 x 1(T 18 


8.9 x 10~ 19 


Haas & Glasssrold (1993) 




1.1 ±0.2 x 10~ 18 




Haas et al. (1995) b 




1.9 ±0.1 x 10~ 18 




ISO Castro-Carrizo et al. (2001) 


[0 I] 145.5 /mi 


11 ±4 x 10- 20 


4.8 x 10- 20 


Haas et al. (1995) 




2.7 ±0.5 x 10~ 20 




Castro-Carrizo et al. (2001) 




5 x 10- 20 




Barlow (1999) (& Priv. Comm) 


[Si II] 34.81 /mi 


0.93 ±0.08 x 10~ 18 


1.5 x 10~ 18 


ISO This paper 




0.94 ±0.37 x 10~ 18 




Haas & Glasseold (1993) 


[C II] 157.7/mi 


1.1 ±0.1 x 10~ 19 


0.7 x 10~ 19 


Barlow (1999) 




1.2 ±0.1 x 10~ 19 




Castro-Carrizo et al. (2001) b 



"Assuming each ion is the dominant ionization state. 
6 Off-source emission reported at 50% level. 
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A. Approximate Absolute Flux Calibration for the TEXES a Ori Spectra 

The analysis of M supergiant atmospheric dynamics can be made without an absolute 
flux calibration of the TEXES spectra, but to explore the full thermodynamic diagnostic 
potential of the [Fe II] 24.52 /im and 17.94 /im, and [Fe i] 24.04 /im emission lines requires 
absolute flux calibrated spectra. This enables a comparison of flux predictions from model 
atmospheres with observations from airborne and space observatories, i.e., [Si II] and [O i] 
emission detected with KAO, and [C II] and [Fe II] emission in ISO spectra. Here we 
describe an approximate flux calibration of the TEXES spectra for a Ori which has several 
independent mid-IR flux measurements that can be used to calibrate and correct for the 
TEXES slit losses. 

In this Appendix we bring the TEXES spectra to an absolute scale by adopting 
the shape of published ISO Short Wavelength Spectrometer (SWS) spectra over the 
16.5-26.5 /zm wavelength region and scale the continuum flux at 25 /im to a value derived 
from a combination of fluxes from color-corrected photometry from the Diffuse Infrared 
. Background Experim ent (DIRBE) on the Cosmic Background Explorer (COBE) satellite 



Boggess et al- 



ii^!), color-corrected InfraRed Astronomic al Satellite (IRAS) 



and the cryogenic grating spectrometer (CGS) on the KAO (IHaas fe 



Glassgold 



are also checked against 8-13 /im UKIRT CGS3 spectrophotometry (IMonnier et al 



photo metry, 



1993 ). Th ese 



19981). 



a Ori's mid-IR spectra contains emission from close to the star and from spatially 
extended optically thin silicate dust emission. The TEXES, ISO, IRAS, KAO-CGS, 
and DIRBE observations all have different entrance apertures and beam sizes and as a 
consequence the TEXES spectra are scaled with differential corrections that account for the 
different TEXES slit losses for the extended dust and point source stellar emission. The 
individual steps are outlined below. 
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A.l. ISO Spectrometers 



To derive color-corrected fluxes for the DIRBE and IRAS photometry the spectral shape 
and system responses are required across each photometric passband. To find the color- 
corrections K\, where the observed color-corrected flux is given by F obs = F DIRBEIR as/K\, 
we adopted ISO spectra and MARCS models (see Table [9]). 



Three different reductions of a Ori's ISO SWS spectra, ob 



resolution o 



(120061), and 



R ~ 1000 a t scan speed #4, from 



Sloan et al 



; ained with a grating 



Justtanont et al. 



(119991 ) 



Verhoelst et al. 



(120031 ) were used to derive K\ for the DIRBE 3.5 fim, 4.9 /im 



12 /im, and 25 /im and IRAS 12, and 25 /im fluxes. T he LWS 



grating spectrurdj which was 



19991 ) was used to find K\ for 



obtained two days before the end of the ISO mission ( IBarlow 
the 60 and 100 /im fluxes. 

The [Fe II] TEXES observations correspond to the ISO SWS band 3 which has a 
14 x 27" aperture. Betelgeuse is a bright IR source and provides an I SO calibration 



chal 



enge in addition to ISO's known systematic flux calibration errors ( Verhoelst et al. 



20061 1 . This is a reason why the published flux spectrum values differ by as much as 40%. In 



the following we boot-strap the ISO flux calibration using the color-corrected photometry. 



A. 2. DIRBE Photometry 



We processed the DIRBE photo metr y from the Calibra ted Individual Observations 
"thl 



in a fashion similar to 



Smith! (12003 ) and 



Smith et al 



(120041 ) . Flux outliers caused by 



cosmic rays coincident with, and off, the source position were rejected. We did not use IR 
databases to reject certain scan directions to avoid potential source confusion which can be 



7 Kindly provided by M. Barlow. 
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important at shorter wavelengths. The co 
are essentially identical to those found by 



or-uncorrec ted flu xes are given in Ta ble [H These 



Smithl 120031) and 



Smith et al 



((20QJ). There are 



typically about 300 measurements in each band. 

At 25 jum and 60 /im the mean of the individual flux error estimates is similar to the la 
standard deviation of the total dataset, while at 12 /im the la standard deviation is three 
times the mean error hinting at intrinsic short term variability during this epoch. Over 
the 10 months of cryogenic DIRBE observations there was a 7% increase in flux at 12 /im 
and a 3% increase at 25 fim. Note that the standard deviation of the mean is an order of 
magnitude smaller than the standard deviation of the individual measurements. 

To convert the DIRBE photometry into color-corrected fluxes the spe ctral distribut ions 



across each band were combined with the spectral response curves from 
(|19981 ) . The color-corrected fluxes are given in Table [9] 



Hauser et al. 



A. 3. IRAS Photometry and Low Resolution Spectrometer 

Betelgeuse had IRAS 12, 25, 60, and 100 /im fluxes measured three times in March 1983 
and the differences in the individual measurements were consistent with their uncertainties. 
The IRAS Point Source Catalog measurements are given in Table O 

The IRAS Low Resolution Spectrometer (LRS) spectra have a resolution of R ~ 20 
and cover 7.5-22.5 p. T he a Ori LRS spectra were combi ned with the correction factors 



from 



Cohen et al. 



IRAS Explanatory Supplement 



(119921 ) and the system responses from the 
(119881 ) to derive a color-corrected 12 /im flux of 3393 ± 136 Jy. For the 25, 60 and 100 /im 
fluxes we used the non-contemporaneous ISO SWS and LWS spectral shape to derive 
color-corrected fluxes which are given in Table [9j 



-67- 



A.4. Kuiper Airborne Observatory 



Haas k, Glassgoldl (119931 ) report the 1992 January 16 detection of [Si II] 34.81 jum and 



[0 i] 63.18 /an emission lines using the cryogenic grating spectrometer (CGS) on-board 
the KAO with resolutions of R ~ 2900 (44" aperture) and R ~ 3700 (34" aperture), 
respectively. The measured line fluxes are: [Si II] = 0.94 ± 0.37 x 10~ 18 W cm -2 and [O i] 
= 2.37 ± 0.21 x 10 -18 W cm -2 . They also measured nearby continuum fluxes of 724 ± 29 Jy 
at 35 fim, and 96 ± 41 Jy at 63 /im. 



A. 5. Adopted Fluxes 



The IR observations discussed abov e were made a t diffe rent epochs and stellar 



variability is an important consideration. 



Monnier et al. 



(119981 ) present seven mean 8-13 /im 



fluxes obtained between 1994 August and 1996 September and the variation of these fluxes 
were within the absolute flux calibration uncertainties of 5-10%. In the mid-1990's there 
is a hint that the 9. 7 /on silicate emi ssion feature shifted slightly to shorter wavelengths 



than in the 1983 IRAS LRS spectra (IMonnier et al. 



19991 ). These authors also found that 



the LRS spectra are somewhat bluer than both previous and later observations, indicating 
a residual miscalibration of the IRAS LRS spectra. In light of this we redetermined the 
color-correction for the IRAS 12 /x m flux (also given in Table [HD using the ISO-SWS spectra 



and this leads to a 15% difference. 



Bester et al. 



(119961 ) give eight 11.15 /im fluxes obtained 



between 1989 November and 1995 August and the variations of these are also consistent 
with a 5% uncertainty. Although the observations are scarce, it appears that intrinsic 
star+dust flux variations at 12 /im on decadal timescales are less than 10%, and at 25 /im 
probably less. 



The DIRBE 25 /im and 60 /im fluxes from 1989-1990 are significantly larger than 



Flux (Jy) 




Fig. 9. — ISO spectra and photometry of a Ori described in Appendix A. The la error 



bars are mostly within the filled circles. The flux curve of Uusttanont et al.l (119991 ) is in 
good agreement with the IRAS fluxes and the DIRBE 4.9 and 12 /im fluxes. The DIRBE 
25 and 60 /xm fluxes are upper-limits because of the large beam size and the presence of 
extended IR emission sur rounding Betelgeuse. The KAO 3a absolute flux un c ertain ty is 25% 



( jHaas fc Glassgold 



Verhoelst et al. 



(120061 ) appears 



19931 ). We note that the calibration of 
systematically low (> 20%) at these wavelengths, but these authors note that their adopted 
multiplicative factors are lower than typically adopted. The mean 8-13 ^im fluxes from 



19m 



UKIRT CGS3 spectrophotometry between 1994 Aug and 1996 Sept (IMonnier et al. 
ranged from 4602-4943 Jy with 1 0% uncertainty in the absolute flux calibration. The mean 



1997 8-13 value derived from 



Justtanont et al. 



(J1999J) is 10% lower at 4386 Jy. A color 



version of this figure is available in the electronic edition. 
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the 1983 IRAS fluxes, and the 25/zm fluxes are greater than the IVerhoelst et al.l (120061 ) 
reduction of the 1997 ISO observations, which is expected to have a typical photometric 



uncertainty of ~ 10% (IVan Malderen et al. 



20041 ) . 



The majority of the stellar wind dust em ission at 25 /im is expected to lie within V 



of a Ori (IStencel. Pesce fc Hagen Bauer 



19881 ). i.e., well within the large DIRBE beam of 



0.°7 x 0.°7 and the IRAS detectors FOV 0.75 x 4.6' . Only a few percent of dust emission is 



Noriega-Crespo et al. 



(119971 ) discovered 



not collected in the ISO band 3C and 3D apertures, 
a wind-ISM bowshock located at 5' from the star, and a nearby linear structure. The 
DIRBE fluxes contain contributions from this low surface intensity structured emission. 
At 60 /im the bowshock alone has ~ 30% of the flux from the star+wind. The DIRBE 
color-corrected 60 ixm flux p± 445 Jy is greater than the combined star+bowshock flux 



(INoriega-Crespo et al. 



19971 ) and while the relative contribution of the extended emission 



at 25 /im is likely to be less than that at 60 /im it probably accounts for some of the excess 
DIRBE 25 /im flux. At 12 /im the DIRBE flux should be dominated by the star and 
its dusty wind, and is consistent with the IRAS flux. So although the DIRBE absolute 
photometry is very good, the large beam size makes it less suitable to normalize the ISO 
spectra for A > 12 /im. 

A comparison of the 4.9 and 12 /im DIRBE, the 12 and 25 /im IRAS, and the 35 /im 



KAO fluxes with the different ISO reductions revea 



Justtanont et al 



lower. Indeed 



1999) but not the 



Verhoelst et al 



Verhoelst et al. 



s a g ood overall agreement with 



20061 ) spectrum which is significantly 



(120061 ) noted that they scaled the sub-band fluxes down 
significantly more than expected for A > 4 /im. Our results suggest that the ISO calibration 
for bright IR sources requires improvement. 

In summary, for the absolute flux scaling of the TEXES spectra we adopt the ISO SWS 
spectrum normalized to the mean of the 12 & 25 /im IRAS fluxes and the 4.9 and 12 /im 
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DIRBE fluxes. We assign an absolute flux uncertainty of 20% to account for the combined 
stellar variability and the scatter in the different mission normalizations. The ISO spectra 
scaled to these fluxes is shown in Figure M 

The relative fluxes for the TEXES [Fe II] lines ultimately relies on the normalization 
of the slightly overlapping ISO spectra in band 3C (A = 16.5 — 19.5 /xm) and band 3D 



(A = 19.5 — 27.3 /mi), e.g., between (A19.37 — 19.57 /mi) (ISloan et al 



2003) 



Verhoelst et al. 



'he difference 



6) 



(120061 ) is 



in the band 3C and band 3D multiplicative factors listed by 
~ 1.3%, we therefore adopt a 2% relative flux error in the continua near the [Fe II] 17.98 //m 
and 24.52 /mi lines. 



A. 6. Putting it all together 



The inner dust radius is measured to be ~ 1" (IDanchi et al.lll994t ) and is comparable 
to the TEXES slit width in the dispersion direction. The observed TEXES spectra thus 
suffer different slit losses for the point source photospheric and chromospheric emission 
and for the more diffuse extended dust continuum emission. To place the TEXES spectra 
onto an absolute flux scale we must first apply corrections for the emission not transmitted 
through the 2 x 17" slit and telescope and instrument losses. To estimate the separate slit 
losses of star and dust we use the silicate dust specific intensity model from HBL01 and 
convolve the resulting sky-image with a Gaussian to represent a combination of seeing, 
diffraction and pointing jitter, whose width is estimated by the recorded TEXES spatial 
profile: FWHM~ 2.6" at 24.5 /mi, and FWHM~ 2.1" at 17.9 /an. 

The flux recorded by TEXES 

FtEXES = Ai oss {C*F* + CdustFdust] (Al) 

where Ai oss is a multiplicative factor to account for combined telescope and instrument light 
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losses not already corrected for in the radiometric flat-field procedure, and C* and Cdust are 
the fractions of the total star and dust flux that pass through the slit, respectively. These 
are calculated assuming the HBL01 sky intensity model, i.e., at 24.5/im: = 0.65 and 
Cdust — 0.31. 

We assume that the ISO aperture records the total flux from the system, i.e., the star 
and wind emission but not the bowshock emission then 

F = F, + F dust . (A2) 

At 24.5 fim, 17.9 /im and 11.15 fim the ratio of dust to star emission derived from the 
HBL01 model is Fdust/F* = 1-6, 1.4,0.7, respectively. Coefficients Ai oss can then be found 
and the TEXES spectra corrected and scaled to the ISO spectrum (which is scaled to the 
DIRBE and IRAS fluxes) using the four TEXES orders recorded at each wavelength setting. 
The resulting TEXES spectra for a Ori are shown in Fig. [21 

We note that the 6" nod of the star along the slit followed by a subtracted image which 
cancels the sky noise clips some of the dust emission which is present at 3". The spatial 
profile shows some emission beyond the Gaussian core and this is also predicted from the 
HBL01 dust model. While this does not affect the present emission line analysis we estimate 
that the flux measurement procedure underestimates the total flux by ~ 5 — 8%. 



B. Composite Model Atmosphere 



To determine the formation radii of Betelgeuse's mid- and far-IR emission lines 
discussed in §5 requires a comprehensive model that encompasses the chromosphere, inner 
wind, and CSE. Curre ntly no such compr ehensive models exist. Models do exist for th e 



inner region, (HBL01: 



Harper et al. 



20011), and the CSE (RG91: 



Rodgers & Glassgold 



19911 ) . and here we describe a composite dynamic and thermodynamic 1-D model that 
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utilized these models, and interpolates between them. 



The HBL01 model was based on the Hipparcos a Or i distance of 131 pc, but fortuitously 



the revised distance of 197 ± 45 pc (IHarper et al. 



20081 ) is also that originally adopted in 



RG91 (200 pc). We therefore take Rodger & Glassgold's stellar pa rameters as our nominal 



Harper et al 



200 



) the Infrared Spatial 



19961 ) was adopted, but it 



values: R* = 1078 R Q , and 0* = 50 mas. In 
Interferometer 11.15 /im angular diameter of 56 mas (iBester et al 
now appears that this may be an over-estimate of the photospheric size (iPerrin et al. 
and the RG91 value is probably closer to the actual value. 
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B.l. Thermal Structure 



The HBL01 model was based on angular resolved radio emission and the old Hipparcos 
distance of 131 pc. Although this model is insensitive to the uncertain angular diameter of 
the photosphere, th e thermodynamic p roperties must be scaled to th e improved distan ce 



estimate of 197 pc (IHarper et al. 



20081 ). The radio interferometry of 



Lim et al. 



(1998) 



essentially measures the angular distribution of the specific intensity. 

If we consider the model in terms of the normalized radial distance Z, where R = ZR*, 
then T gas (Z) remains unchanged but .R* is now a factor of 197/131 larger. The angular 
radio brightness distribution also requires that the optical depth through the atmosphere 
remains unchanged, so that r oc n e nndR is constant. Therefore to satisfy this constraint 
the particle densities each scale as 1/ a/197/131. 

The single component 1-D temperature structure derived from the radio represents 
a complicated averaging of the electron temperatures of the hot chromospheric plasma 
and cool wind plasma. We expect that the filling factor of the hot plasma decreases with 
increasing radius, so the bulk of the plasma is cooler than inferred from the radio. For the 
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calculation of the mid-IR emission we adopt a lower temperature distribution that joins the 
HBL01 model at Z = 7. 

When the RG91 models were constructed it was widely believed, on the basis of 
theoretical grounds and semi-empirical models based on spatially unresolved data, that 
the inner wind had warm chromospheric temperatures and RG91 adopted a nominal inner 
boundary condition (BC) of T gas = 8000 K at Z = 3. They also provided variational 
calculations for the temperature structures resulting from different mass-loss rates and 
where the inner BC was set to T gas = 4000 K. In contrast T gas = 2764 K for HBL01. The 
RG91 models with different inner temperature BC's have similar shapes and smoothly 
converge to join at 30-R*. Interior to Z = 7 we take the temperature structure from HBL01 
which is constrained by the long-wavelength radio observations, and exterior to that the 
temperature structure is obtained by extrapolating on the difference between the RG91 
models with inner BC of T e = 4000 K and 10,000 K. 

The CSE temperature structures are less sensitive to differences in mass-loss rates: the 
RG91 value of 3.5 x 10~ 6 MQyr^ 1 (assuming a mean mass per hydrogen nuclei S = 1.4 
with V W i n d = 10 kins -1 ) is similar to the rescaled value of HBL01 (see below). The 
adopted composite temperature structures are shown in Figure [10J The solid line includes 
contribution from hot plasma while the dashed line is a schematic representation of the 
temperature of the cool wind. 

B.2. Velocity Fields 

A detailed description of the run of the mean outflow velocity, V(R), and fluctuations 
about this value, Vt ur b{R), are important for both emission line profile calculations and for 
determining the photon escape probability which enters into the flux emitted by the envelope. 




Fig. 10. — Composite temperature structure for a Ori's chromosphere, inner wind and CSE. 



The inner region (solid line) is based on a distance- scaled HBL01 spatial l y-exte nded semi- 

Rodgers fc Glassgo J ( 1991 ). The dash 



empirical model, and the CSE is based on a model of ; 
line is a schematic representation of the temperature of the dominant plasma (see text for 
details.) 
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Vturb(R) which enters into the line profile calculations is assumed to be random and isotropic 
(in the absence of any other knowledge) and given by V turb (R) = ^V t 2 herm + V^ on _ therm , i.e., 
it includes the thermal motion of the diagnostic species combined with, as yet unidentified, 
non-thermal motions. Here V non - t herm dominates. 



B.2.1. Mean Radial Outflow Velocity 



Goldberg! (119791 ) discussed the empirical wind velocity constraints from atomic and 
molecular absorption lines. There are two clearly identified velocity features and a dynamic 
region which may be physically distinct. The S2 shell which is narrow and discrete in 
velocity-space [Vt ur b < 1 km s _1 and V(S2 ) g 17 km s -1 ] is observed in both absorption 



and scattered emiss ion from atomic lin es (IMauron 



fundamental band (IBernat et al 



19901 ). in absorption in the CO 4.6 /xm 



19791 ). and in emission in millimeter CO rotati onal 



emission. In K I the S2 shel 



inner edge of ~ 7" (IMauron 



is seen out to a radius of 50" (IPlez &: Lambert 



20021 ) with an 



19901 ). However, the similarity of the shape of CO mm-radio 



emission profile in single-dish observations w i th di fferent beam-sizes suggests the outer 



edge of the CO is < 12", e.g., 



inner edge of 2.5" is plausible (IHugg 



Hug gins et al 
insf" 



(1994J ), and from modeling considerations an 



19871 ). The spatial extent of different diagnostics 



is unlikely to be identical because of differing ionization balances in the extended envelope. 
We shall find that the atomic and singly ionized CSE emission lines are formed within 
these inner radii and are likely to be characteristic of the SI shell [Vt ur b ~ 4kms~ 1 and 
V(S1) ~ 10 km s -1 ]. SI is observed in P-Cygni profiles and is blended with the photospheric 
absorption line, which may lead to a small over-estimate of the radial velocity of the shell. 
While SI has not yet been spa tially resolved in m m-radio CO observations, which show 



that it lies within the S2 shell (IHarper et al. 



20091 ) . t he SI shell has bee n resolved in the 



photospheric scattered CO 4.6 /xm fundamental band (ISmith et al. 



2009|). 
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Dynamic flow features have been observed in ultraviolet (UV) Fe II line profiles 



(IBoesgaard fe Magnan 



1975|), and 



Carpenter! (119841 ) schematically mapped out the radial 



wind velocity using observations from the International Ultraviolet Explorer. The wind 
acceleration is mo re seen clearly in the G o ddard High Resolution Spectrograph (GHRS) 



spectra studied by 



Carpenter fe Robinson! (119971 ) (see their Figure 6). The line profiles map 



the flows close to the star and this acceleration region may be part of SI s tructure. At 



1-1.5" CO 4.6 Mm wind scattering observations from Phoenix on Gemini-S (IHarper et al. 



2009 



Smith et al. 



20091 ) reveal that the wind has the velocity of the SI shell. The emission 



for most CSE lines (except carbon) is formed interior to 1" so we will adopt a wind model 
that reaches a terminal speed of ~ 10 km s" 1 at 1" and extended out beyond 5". Given that 
it is unlikely that there is a smooth transition in the wind properties between the SI and 
S2 material this is a reasonable procedure, until further spatial information is obtained in 
this interesting region. 

We estimate the run of wind velocity close to the star using the new temperature 
structure to locate the normalized radius (Z) where the absorption minima (V a b s ) of UV 
GHRS Fe II wind features have a radial optical depth r = 1. We then select a wind velocity 
profile that approximates a range of Fe II Z — V^, values and that is also the solution to 



the constant pressure wind equation (IBrandt 



197(1 ) 



V 2 

^crit 



In 



V 2 



R Rcrit 



R 



crit 



R 



3. 



(Bl) 



The velocity profile is the solution of this transcendental equation, and is defined by 
V cr it = 2.5 km s _1 and Z cr i t = 2.75. We also limit the wind speed not to exceed 10 km s _1 at 
large distance. 



-77- 



B.2.2. Turbulence Velocities 



For the radial distribution of turbulence, we note that the studies of eclipsing binaries 



reveal that the turbulen t velocities in the chromosphere and inner wi nd are typica 



hydrogen sound speed (jEaton 



1993 



Baadeetal 



19961 ). but see also 



Kirsch et al- 



l y ~ 1. 5 x 



(2001) 



These small scale mot ions are proba bly related to the mass-loss mechanisms and may reflect 
MHD waves, e.g., see 



Jordan! (119861 ). Here we adopt ~ 1.5 X hydrogen sound speed, namely 

Vnon-therm (R) = 0.19^T gas (R) km S" 1 . (B2) 

We note that this gives V tur b — Vnon-therm ~ 12 km s -1 near the base of the wind in 
agreement with a Ori's [Fe II] line widths. In this model at Z ~ 9 the wind and non-thermal 
turbulence velocities are approximately equal. 



B.3. Density Structure 

For the hydrogen densities we join the distance-scaled HBL01 model with the RG91 
model. The densities in the RG91 model are oc M/V^ and they originally adopted a 
terminal wind speed of V^, = 16 km s" 1 (i.e., the S2 velocity), however, it now appears 
that the appropriate win d speed for the most of the CSE line formation is 10 km s _1 



( iHaas &: Glassgoldl Il993l ). To maintain the same density structure with this lower 
velocity requires lowering the original RG91 mass- loss rate of M = 5.6 x 10" 6 M^r" 1 to 
M = 3.5 x 10 -6 M Q yr _1 . We define the hydrogen density for Z > 7 assuming 

M 

HH ~ 4n&m H i;V(R) ' (B3) 
The densities are therefore increased over the constant wind velocity limit. At Z = 7 the 
inner densities are a factor ~ 2 larger than implied by Eq IB3I so the two have a simple 
join. The mass-loss rate implied by the UV Fe II lines with this new density structure is 
4.8 ± 1.3 x lO _6 M yr _1 which is consistent with the lower velocity RG91 value. 
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B.4. Inhomogeneities 



The 1-D thermal structure derived from radio interferometry represents a mean value of 
the different structures that co-exist at a given stellar radius. Near the temperature peak a 
crude estimate gives an area filling factor of hot chromospheric material of A r h rn m. (Z) ~ 1/ 3 



Harper et al. 



20011 ). while by Z ~ 3 the filling factor is much smaller (lHarper fc Brown 
20061 ). The filling factor of the hot plasma is expected to continue to decrease with 
increasing radius. The temperature of the bulk plasma in the region that encompasses the 
hot chromosphere might then be lower than the HB01 model, and an alternate schematic 
model is shown as a dash line in Fig. [TO], and this is adopted in §5 to calculate the 
contribution functions. 
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Table 9: DIRBE and IRAS photometry: color-uncorrected and color-corrected flxues. 
Band Wavelength Color-Uncorrected Color-Corrected K\ Spectral shape a 



Flux ±lci (Jy) Flux ±1(7 (Jy) 



DIRBE 














Band 1 


1.25 


3.04 ±0.11 x 


10 4 


3.20 ±0.12 x 


10 4 


0.95 MARCS 


Band 2 


2.2 


3.03 ±0.09 x 


10 4 


3.62 ±0.10 x 


10 4 


0.88 MARCS 


Band 3 


3.5 


1.71 ±0.05 x 


10 4 


1.78 ±0.05 x 


10 4 


0.96 ISO-SWS 


Band 4 


4.9 


7.16 ±0.31 x 


10 3 


7.46 ±0.32 x 


10 3 


0.96 ISO-SWS 


Band 5 


12 


3.86 ±0.14 x 


10 3 


3.68 ±0.16 x 


10 3 


1.03 ISO-SWS 


Band 6 


25 


2.91 ±0.06 x 


10 3 


1.48 ±0.04 x 


10 3 


1.95 ISO-SWS 


Band 7 


60 


6.29 ± 1.02 x 


10 2 


4.31 ±0.70 x 


10 2 


1.46 ISO-LWS 


Band 8 


100 


1.02 ±2.76 x 


10 2 


0.87 ±2.36 x 


10 2 


1.17 ISO-LWS 



IRAS 



12 

25 
60 
100 



4.68 ± 0.19 x 10 2 



3.39 ± 0.14 x 10 3 1.38 IRAS LRS 
3.87 ± 0.16 x 10 3 1.21 ISO-SWS 



1.74 ± 0.07 x 10 3 1.20 ± 0.05 x 10 3 1.45 ISO-SWS 
2.99 ± 0.21 x 10 2 2.23 ± 0.16 x 10 2 1.34 ISO-LWS 
9.59 ± 1.92 x 10 1 9.13 ± 1.83 x 10 1 1.05 ISO-LWS 



a The spectral shapes adopted to find the color-corrections, K\. 



